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RESEARCH  SUMMARY 

The  work  performed  under  this  contract  has  had  a  significant 
impact  on  existing  abilities  to  predict  the  behavior  of  tropical  cyclones 
and  on  abilities  to  model  the  predicted  behavior  of  high  altitude  nuclear 
explosions  (HANEs).  This  work  has  resulted  in  the  generation  of  four 
pertinent  publications,  each  of  which  is  described  below  and  attached 
as  an  appendix. 

In  "The  Impact  of  Satellite-Sensed  Winds  on  Intensity  Forecasts 
of  Tropical  Cyclones",  it  is  indicated  that  assimilation  of  these  winds 
with  marine  surface  (or  low-level)  wind  alone  does  not  improve  intensity 
forecasts  appreciably.  A  strong  relaxation  coefficient  in  the  initiali¬ 
zation  scheme  causes  model  rejection  of  the  assimilation  and  an  attenuating 
relaxation  coefficient  is  recommended.  However,  when  wind  observations  at 
the  outflow  level  are  included  in  the  assimilation,  forecasts  improve 
substantially  and  the  best  forecasts  are  achieved  when  observations  over 
the  entire  lower  troposphere  are  assimilated. 

In  "Numerical  Simulation  of  the  Influence  of  Sea-Surface  Temperature 
on  Translating  Tropical  Cyclones",  it  is  suggested  that  the  sea-surface 
temperature  distribution  not  only  affects  the  intensity  and  path  of  tropical 
cyclones  frictional ly,  but  also  thermally.  The  enhanced  evaporation  and 
convergence  over  the  warm  SST  provide  a  favorable  condition  for  the  growth 
of  the  tropical  cyclone,  and  leads  to  a  gradual  shift  of  the  storm  center 
toward  the  warm  ocean. 

"Test  of  a  Planetary  Boundary  Layer  Parameterization  Based  on  a 
Generalized  Similarity  Theory  in  Tropical  Cyclone  Models",  shows  that  in 
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C  comparison  with  a  complex,  one-dimensional,  multi-layer  PBL  model,  the 

GST  parameterization  yields  accurate  moisture  fluxes,  but  slightly  over¬ 
estimates  the  momentum  flux  and  underestimates  the  sensible  heat  flux, 
t  The  GST  parameterization  produces  very  realistic  dynamics,  energetics 

and  thermal  structure  in  an  axisymmetric  tropical  cyclone  model. 

In  the  HANE  phenomenology  portion  of  this  research  effort, 

*  considerable  effort  was  devoted  to  enhancing  and  improving  the  capability 

of  the  PHARO  simulation  code.  These  improvements,  as  code  modifications, 
have  been  delivered  to  NRL  and  are  now  integrated  into  PHARO.  Another 
portion  of  this  effort  led  to  a  unique  interpretation  of  photographic 
data  acquired  on  the  STARFISH  high  altitude  nuclear  event  and  this  work 
is  being  published  separately  as  an  NRL  Memo  Report  because  it  is  classified. 
The  title  of  this  report  will  be  “Interpretation  of  High  Altitude  Streaks 
in  Starfish  (U),  authored  by  S.  Brecht,  C.  Agritellis,  E.  Hyman  and 
P.  Palmadesso.  The  report  will  be  published  at  the  Confidential  level. 
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THE  IMPACT  OF  SATELLITE-SENSED  WINDS 
ON 

INTENSITY  FORECASTS  OF  TROPICAL  CYCLONES 
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ABSTRACT 

The  impict  of  the  satellite-sensed  winds  on  the  intensity  forecasts 
of  tropical  cyclones  is  evaluated  by  a  simulation  study  with  an  axisym- 
metric  numerical  model.  The  parameterized  physics  in  the  forecast  model 
are  deliberately  made  different  from  those  in  the  model  that  generates  the 
observation.  Model  generated  "observations"  are  assimilated  into  fore¬ 
casts  by  12  h  dynamic  initialization. 

A  series  of  24  h  forecasts  with  and  without  assimilation  of  satellite- 
sensed  winds  are  conducted  and  compared  with  the  observations.  Results 
indicate  that  assimilation  with  marine  surface  (or  low-level)  wind  alone 
does  not  improve  intensity  forecasts  appreciably,  that  a  strong  relaxation 
coefficient  in  the  initialization  scheme  causes  model  rejection  of  the 
assimilation,  and  that  an  attenuating  relaxation  coefficient  is  recommended. 
However,  when  wind  observations  at  the  outflow  level  are  included  in  the 
assimilation,  forecasts  improve  substantially.  The  best  forecasts  are 
achieved  when  observations  over  the  entire  lower  troposphere  are  assimilated. 

Additional  experiments  indicate  the  errors  in  the  satellite  observations 
contaminate  the  forecast.  But  the  assimilation  of  inflow  and  outflow  winds 
still  improve  the  intensity  forecast  if  the  satellite  observation  errors 
are  less  than  or  about  the  same  magnitude  of  those  in  the  initial  wind 
field. 
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Introduction 


The  improvement  of  forecast  skill  in  tropical  cyclones  evident  in 
the  1960's  has  not  been  continued  in  the  1970's  in  spite  of  improved  tech¬ 
nology  and  continuing  effort.  The  lack  of  improvement  has  been  attri¬ 
buted  to  the  imperfect  knowledge  of  the  initial  fields  for  objective  models. 
Elsberry  (1977)  attributed  the  poor  performance  of  his  prediction  model 
to  the  deficiency  of  the  initial  wind  data.  For  the  1976  Atlantic  tropical 
cyclone  season,  Hovermale  and  Livezey  (1977)  showed  the  errors  for  the 

36-  and  48-h  track  forecasts  i  lcreased  by  approximately  a  factor  of  three  for 
storms  over  data-void  ocean  regions  as  compared  to  storms  near  coastal 
stations.  In  addition,  the  theory  of  geostrophic  adjustment  requires 
that  the  mass  field  adjusts  to  the  momentum  field  for  low  latitudes  and 
tropical  cyclone  scales  of  motion . (Monin  and  Obukhov,  1959;  Washington, 
1964).  It  is  beyond  doubt  that  wind  observations  are  essential  to  tropical 
cyclone  forecasts. 

Although  initial  wind  analyses  of  tropical  cyclones  have  been 
improved  from  reconnaissance  aircraft  flights,  the  quantity  and  especially 
the  aerial  coverage  of  wind  data  so  obtained  are  inadequate  for  numerical 
model  initializations.  Remote  measurements  from  geostationary  and  orbiting 
satellites  will  be  relied  upon  as  important  data  sources.  Rodgers  et  al 
(1979)  explained  techniques  to  derive  low-level  and  outflow  level  winds 
for  tropical  cyclones  by  tracking  clouds  using  successive  satellite  images. 
Their  results  are  encouraging  in  spite  of  some  difficulties  such  as  the 
short  lives  of  cloud  turrets  and  overcast  conditions  near  the  storm  centers. 

An  experimental  oceanographic  satellite  known  as  SEASAT-1  during  its  short 


lifetime  provided  an  additional  data  source  to  define  large  and  mesoscale 
wind  fields  near  tropical  cyclones.  A  special  scatterometer  (SASS)  flown 
with  SEASAT-1  measured  the  marine  surface  microstructures  and,  through 
appropriate  algorithms,  marine  surface  winds  were  inferred. 

The  purpose  of  this  study  is  to  evaluate  the  impact  of  the  satellite- 
sensed  winds  on  tropical  cyclone  intensity  forecases.  As  a  simulation  study, 
data  generated  by  numerical  models  will  be  used  in  place  of  real  data.  The 

general  strategy  of  such  simulation  studies  follows  that  of  Charney  et  al 
(1969).  First,  a  control  integration  of  the  numerical  forecast  model  is 
performed  to  generate  the  "true"  history  of  the  atmosphere  or  the  "observa¬ 
tion".  A  series  of  "standard  forecasts"  is  then  generated  based  on 
different  initial  states.  Finally,  a  series  of  forecasts  with  the 
"observations”  assimilated  are  conducted  to  evaluate  the  impact  of  assimila¬ 
tion.  For  detailed  reviews  of  such  simulation  studies  and  their  general 
strategies,  readers  are  referred  to  McPherson  (1975)  and  Bengtsson  (1975). 

As  a  preliminary  study,  an  axisymmetric  tropical  cyclone  model  is 
employed  in  this  study,  therefore  only  the  impact  on  intensity  can  be 
studied.  The  method  of  assimilation  used  is  the  dynamic  initialization 
by  relaxation  (DIR)  technique.  To  approximate  the  relationship  between 
the  real  atmosphere  and  forecast  models,  parameterized  physics  in  the 
model  that  generates  forecasts  (forecast  model)  are  altered  from  those 
in  the  model  that  generate  the  observations  (natural  model).  In  the 
following  sections,  the  numerical  model,  the  experimental  design,  and  the 
method  of  assimilation  will  be  discussed  in  sequence.  Finally,  the 
results,  conclusions  and  proposed  future  research  will  be  presented. 
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2.  Numerical  Model 

The  axi symmetric  tropical  cyclone  model  used  in  this  study  is 
similar  to  the  one  described  in  Chang  (1977)  and  Anthes  and  Chang  (1978). 

The  governing  equations  are  in  primitive  form  and  are  in  o(=P/P$)  coordinates 
The  explicit  water  vapor  cycle  and  parameterization  of  cumulus  convection 
follows  Kuo  (1974)  and  Anthes  (1977).  The  boundary  layer  is  contained  in 
the  lowest  model  layer  and  parameterization  of  various  vertical  fluxes  is 
ba:,ed  on  a  generalized  similarity  theory  in  which  Yamada's  (1975)  universal 
functions  are  used  (Chang,  1980).  Charnock's  equation  is  applied  to  com¬ 
pute  marine  surface  roughness  length. 

The  model  atmosphere  is  divided  into  six  layers  (Table  1).  A  uni¬ 
form  horizontal  grid  interval  of  30  km  is  used  from  the  center  to  a  radius 
of  600  km.  The  grid  interval  is  progressively  increased  by  a  factor  of 
two  outside  600  km.  The  leapfrog  temporal  integration  method  with  the 
time-averaged  pressure  gradient  force  (Brown  and  Campana,  1978)  is  employed 
for  numerical  integration.  The  spatial  finite  differencing  is  of  the  second 
order.  The  mean  hurricane  season  sounding  (Sheets,  1969)  is  used  for  the 
initial  and  lateral  boundary  conditions.  The  coriolis  parameter  f  has  the 
constant  value  5  x  10  s  . 

A  period  of  36  h  of  the  control  run  during  which  the  model  tropical 
cyclone  undergoes  a  rapid  intensification  is  chosen  as  the  "truth"  or 
"observation"  (thereafter  referred  to  as  such  in  this  study  except  as  stated 
otherwise).  For  convenience,  -12  h  and  24  h  are  designated  as  the  start 
and  end  of  this  period.  We  selected  a  period  of  rapid  intensification  for 
study  in  order  to  magnify  errors  in  the  forecasts. 
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3.  Experimental  Design 

In  recent  years,  many  simulation  studies  have  been  conducted  to 

evaluate  the  impact  of  incomplete  observation  data  on  numerical  predictions 
(e.g.,  Charney  et  al_,  1969;  Kasahara  and  Williamson,  1972;  Morel  and 
Talagrand,  1974;  Anthes,  1974;  Cane  et  a^,  1979).  In  a  similar  manner, 
numerical  integrations  conducted  for  this  study  can  be  grouped  into  three 
components  (Table  2): 

(1)  Nature  run  -  A  !6  h  segment  of  life  history  of  tropical 
cyclones  designated  as  observation  as  defined  in  Section  2. 

(2)  Standard  foreca; ts  -  A  36  h  forecast  starting  from  -  12  h 
and  a  24  h  forecast  starting  from  0  h.  The  initial  conditions 
for  standard  forecasts  are  obtained  by  the  static,  non- 
divergent  initialization  method  based  on  the  nature  run. 

(3)  Forecasts  with  assimilation  -  12  h  preforecast  integrations 
starting  from  -12  h,  during  which  satellite-sensed  winds 
are  assimilated  into  the  model  solution  followed  by  24  h 
forecasts  starting  from  0  h. 


A  unique  characteristic  of  previous  simulation  studies  is  that  a 
prediction  model  will  make  an  error-free  forecast  given  error-free  model - 
generated  initial  conditions  (.Williamson,  1973).  Because  this  is  rather 
unrealistic,  errors  of  various  kinds  were  added  to  the  observations  either 
in  initial  conditions  for  the  forecasts,  or  in  data  for  assimilation. 

Both  random  errors  (e.g.,  Williamson  and  Kasahara,  1971)  and  bias  errors 
(e.g.,  Anthes,  1974)  have  been  introduced  into  the  observations  in  previous 
studies.  Forecast  runs  with  initial  random  errors  sometimes  exhibit 
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Table  2.  List  of  Simulation  Experiments 


v  are  radial  and  tangential  winds,  respectively.  Subscripts  denote  layers. 


unrealistic  error  growth  characteristics  because  gravity  waves  and  model 
physics  act  to  smooth  them.  Besides,  random  observational  errors  are  not 
the  major  problem  with  real  data,  where  systematic  errors  are  known  to 
have  caused  more  problems  (McPherson,  1975).  Biased  errors  are  generally 
determined  subjectively  and  may  be  unwarranted  and  unrealistic.  No  error 
is  artifically  added  in  the  initial  fields  for  all  forecast  runs  in  our 
study,  instead,  errors  in  the  initial  wind  fields  are  introduced  by  the  non- 
divergent,  gradient-balanced,  static  initialization  procedure  adopted  here. 
Such  initialization  procedure  is  currently  in  use  operationally.  Figure  1 
shows  the  errors  of  speed  in  the  initial  wind  fields  of  forecasts  at  -12 
h.  As  expected,  large  errors  occur  in  the  low-level  and  the  outflow  level, 
where  divergent  components  of  wind  vectors  are  largest.  The  initial  errors 
for  forecasts  initialized  at  0  h  have  the  same  characteristics. 

In  previous  simulation  studies,  the  models  which  generated  the  fore¬ 
casts  were  identical  with  the  models  that  generated  the  observations  (see 
McPherson,  1975,  for  review).  This,  of  course,  is  very  unrealistic.  In 
reality,  numerical  forecasting  models  with  finite  spatial  resolutions  and 
parameterized  physics  cannot  reproduce  the  atmosphere  even  if  perfect  initial 
conditions  are  obtained.  To  properly  account  for  the  discrepancies  between 
real  forecasting  models  and  the  atmosphere,  the  parameterized  physics  in 
our  forecast  model  are  deliberately  altered.  The  parameters  changed  are 
those  we  feel  most  uncertain  about  in  the  physical  parameterizations  in  current 
numerical  models,  namely,  the  effective  air-sea  exchange  coefficients  and 
the  vertical  distribution  of  latent  heat.  The  effective  coefficients  of 


eddy  transfers  of  momentum,  sensible  heat  and  latent  heat  in  the  forecast 
model  are  set  at  90X  of  those  in  the  nature  model.  The  lot  error  is  well 
within  the  expected  error  in  the  boundary  layer  formulations.  The  vertical 
distribution  of  cumulus  heating  is  also  changed  so  that  approximately  5%  of 
the  heating  in  the  lower  troposphere  is  shifted  to  upper  troposphere.  The 
5%  difference  is  within  the  variation  of  the  observed  heating  distributions. 
Due  to  these  two  changes,  the  forecasts  in  our  study  do  not  asymptotically 
approach  the  observations  even  after  long  integration.  Because  the  errors  in 
the  parameterized  physics  in  our  forecast  model  are  within  the  differences 
between  the  atmosphere  and  the  current  operational  forecast  models,  the 
asymptotic  difference. between  our  "forecasts"  and  "observation"  is  quite 
realistic. 
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4.  Method  of  Assimilation 

The  satellite-sensed  data  are  in  many  occasions  incomplete  in  that 
they  do  not  contain  observations  of  all  meteorological  variables  simul¬ 
taneously  or  the  observations  are  made  at  different  locations  and  times. 

To  incorporate  such  data  in  a  dynamically  consistent  way  into  the  numerical 
models,  suitable  methods  of  assimilation  must  be  used.  From  the  direct 
insertion  method  (e.g.,  Charney  et  al_,  1969)  to  the  complicated  variational 
assimilation  (e.g.,  Sasaki,  1969),  there  are  many  methods  of  cssimilation 
in  existence.  However,  not  all  methods  are  applicable  for  the  satellite- 
sensed  data  in  question.  The  wind  fields  derived  from  GOES  images  are 
basically  restricted  to  low  and  outflow-levels  in  tropical  cyclones 
(Rodgers  et  al_,  1979),  and  wind  fields  measured  from  SEASAT-1  are  at 
aneomometer  level.  For  such  data  with  poor  vertical  resolution,  a  method 
called  dynamic  initialization  by  relaxation  (DIR)  is  desirable. 

DIR  is  a  technique  wherein  the  meteorological  variables  are  relaxed 
(or  nudged)  by  using  the  model's  governing  equations  toward  the  observed 
values  during  a  preforecast  integration  (Anthes,  1974;  Hoke  and  Anthes, 

1976).  The  technique  has  shown  great  promise  in  real  data  applications 
(Nitta  and  Hovermale,  1969;  Davies  and  Turner,  1977;  Hoke  and  Anthes,  1977). 
Mathematically,  governing  equations  during  the  preforecast  integration  are 
modified  to: 

N 

If*  F  (X,t)  +  £  X(c  ,  fit,  fir,  5z)  (x°  -  x)  (1) 

dZ  -  n=l  n 

where  x  is  an  element  in  the  vector  of  variables  X,  the  function  F  contains  the 
normal  terms  in  governing  equations,  x°  is  the  observation,  N  the  number 
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of  observations,  andAthe  relaxation  coefficient.  In  a  full  four  dimen¬ 
sional  assimilation,  A  is  the  function  of  observational  error,  en,  the 
time  separation  of  the  observation,  fit,  the  horizontal,  (fir),  and  the 
vertical,  (5z),  spatial  separations  between  observations  and  grid  points. 

It  should  also  depend  on  the  meteorological  variables. 

To  simplify  the  functional  form  of  A,  we  will  use  point-to-point 
relaxation,  i.e.,  variables  are  relaxed  toward  observations  made  at  the 
same  model  grid  points  only.  This  requires  that  observations  be  taken  at 
model  grid  points  and  all  observations  be  taken  simultaneously.  Note  that 
the  horizontal  resolution  of  the  satellite  measurements  does  approach  those 
of  typical  operational  forecast  model  of  tropical  cyclones.  With  the  develop¬ 
ment  of  suitable  boundary  layer  models  for  vertical  extrapolation  (Yu,  1980), 
the  convenience  of  point-to-point  relaxation  assumed  for  convenience  in  this 
study  is  nearly  available  in  operational  forecasting.  The  time  lag  of 
measurements  over  the  domain  of  tropical  cyclones  within  one  satellite 
revolution  is  negligibly  short  as  compared  to  the  12-h  preforecast 
integration.  We  take  note  that  the  swath  width  of  orbiting  satellites 
nevertheless  may  not  be  large  enough  to  cover  the  entire  tropical 
cyclone. 

The  satellite-sensed  winds  are  not  free  of  errors.  Rodgers  et  al 
(1979)  estimated  the  mean  speed  errors  in  their  derived  winds  to  be 
2.5  m  s*1  relative  to  aircraft  measurements.  There  are  conflicting 
reports  on  the  errors  of  SEASAT  measurements  CBlack,  1979;  Jones  and 
Pierson,  1978),  but,  in  general,  the  errors  of  satellite-sensed  winds  are 
smaller  than  those  introduced  by  the  objective  analyses  over  the  oceans 
(Cardone  et  al_,  1976).  The  contribution  of  satellite  measurements  is  not 
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in  the  general  error  reduction  but  in  the  filling  of  data-void  areas 
(Ghil  et  al_,  1979).  For  a  clear  demonstration  of  the  impact  in  assimilat¬ 
ing  winds  at  different  levels,  in  Exps.  2-9  it  is  justifiable  to  assume 
that  the  satellite  observations  are  error-free  in  comparison  to  the  initial 
and  model  errors.  However,  errors  of  different  magnitudes  are  added  to 
the  satellite-sensed  wind  in  Exps.  9E  and  9E2  to  evaluate  the  extent  to 
which  the  observation  errors  contaminate  the  forecast. 

The  equations  of  motion  in  the  preforecast  integration  in  this  study 
can  simply  be  written 

It  Vk  =  F  (X,t)  +  A(v°  -  vk)  .  (2) 

where  k  denotes  the  layer  in  the  model  where  observations  are  available. 
Three  different  values  for  A  are  tested:  A^  -  10~4  s_1  for  weak  relaxation, 
A.  =  10~3  s"1  for  strong  relaxation;  and  Aa  =  A„  (6t  -  t)/12,  -12  h<  6t£  0 
for  attenuating  relaxation.  Figure  2  illustrates  the  time  variations  of  A. 
Observations  are  assumed  to  be  taken  at  0  h  and  are  assimilated  into  model 
prediction  during  -12  to  0  h  in  all  of  the  assimilation  experiments 
(Exps.  4-9). 
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in  the  general  error  reduction  but  in  the  filling  of  data-void  areas 
(Ghi  1  et  al_,  1979).  For  a  clear  demonstration  of  the  impact  in  assimilat¬ 
ing  winds  at  different  levels,  in  Exps.  2-9  it  is  justifiable  to  assume 
that  the  satellite  observations  are  error-free  in  comparison  to  the  initial 
and  model  errors.  However,  errors  of  different  magnitudes  are  added  to 
the  satellite-sensed  wind  in  Exps.  9E  and  9E2  to  evaluate  the  extent  to 
which  the  observation  errors  contaminate  the  forecast. 

The  equations  of  motion  in  the  preforecast  integration  in  this  study 
can  simply  be  written 

ft- h  =  F  +  x(y£  -  yk)  .  (2) 

where  k  denotes  the  layer  in  the  model  where  observations  are  available. 
Three  different  values  for  X  are  tested:  X^  =  10~^  s"1  for  weak  relaxation 
Ac  =  10"^  s"1  for  strong  relaxation;  and  X.  =  Xc  ( 6t  -  t)/12,  -12  h<  6t£  0 
for  attenuating  relaxation.  Figure  2  illustrates  the  time  variations  of  X. 
Observations  are  assumed  to  be  taken  at  0  h  and  are  assimilated  into  model 
prediction  during  -12  to  0  h  in  all  of  the  assimilation  experiments 


5. 


Standard  Forecasts 


During  the  period  between  -12  h  and  24  h,  observation  shows  a  rapid 
intensification  of  the  tropical  cyclone,  the  minimum  central  pressure 
deepens  from  998  mb  to  953  mb  (Figure  3)  and  the  maximum  wind  speed  increases 
from  29  m  $~*  to  52  m  s”*  (Figure  4).  A  36  h  forecast  starting  from  -12  h 
(Exp.  2)  and  a  24  h  forecast  starting  from  0  h  (Exp.  3)  are  conducted  based 
on  the  error-free  initial  mass  field  and  non-di vergent,  gradient-balanced 
wind  fields.  Typical  wind  speed  errors  in  the  initial  wind  field  are 
illustrated  in  Figure  1. 

As  expected  from  previous  experience,  both  forecasts  have  an  initial 
dissipation  stage  due  to  the  onset  of  the  surface  friction.  The  weakening 
of  storm  intensity  is  especially  pronounced  in  Exp.  3  in  that  it  has 
larger  intensity  errors  before  12  h  than  Exp.  2  which  is  initialized  12  h 
earlier.  This  is  indicative  of  the  inadequacy  of  the  static,  non-di vergent 
initialization  employed.  Improvement  of  forecast  during  the  initial  hours 
can  be  achieved  by  using  a  divergent  static  (Jarbell,  1979)  or  a  dynamical 
(Hovermale  and  Li vezey^ 1977;  Kurihara  and  Bender,  1979)  initialization 
scheme. 

After  the  radial  circulations  develop,  both  forecasts  reproduce  the 
observed  intensification  but  at  slower  rates.  After  12  h,  the  24  h  fore¬ 
cast  (Exp.  3)  yields  better  prediction  than  the  36-h  forecast  by  approximately 
2  mb  in  minimum  pressure  and  2  m  s~*  in  maximum  wind  speed.  Both  forecasts 
predict  weaker  storm  intensity  as  compared  to  the  observation.  The 
difference  between  the  observed  and  predicted  intensities  at  24  h  is  about 
15  mb  in  minimum  pressure  and  10  m  s"*  in  maximum  wind  speed.  The  divergence 
of  the  forecasts  from  the  observation  is  a  consequence  of  the  "imperfect” 
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physical  parameter! zations  in  the  forecast  model. 

We  selected  the  root-mea-square  errors  (e)  as  a  measurement  of  the 
accuracy  of  the  predictions  (Panofsky  and  Brier,  1968).  Evolutions  of  e 
for  wind  speed  (.V),  tenperature  CO,  and  specific  hunidity  (q)  with  respect 
to  the  observation  for  Exps.  2  and  3  are  shown  in  figures  5,  6,  and  7, 
respectively.  The  initial  e(.V)  is  large  at  about  4ms1  in  both  experi¬ 
ments  due  to  the  non-divergent  initialization.  It  decreases  for  the  first 
18  h  in  both  experiments  as  the  forecast  storms  intensify.  As  evident  in 
Exp.  2  after  6  h,  values  of  e(.V)  begins  to  increase,  indicating  a  deteriorat 
ing  forecast. 

Both  e(T)  and  e(q)  increase  rapidly  after  initialization  with  time  from 
the  error-free  mass  field.  Their  values  escalate  to  1.5°K  in  temperature 
and  0.9—1, 1  g  kg'1  in  specific  humidity  at  24  h.  Here  again,  Exp.  3 
produces  a  better  prediction  than  Exp.  2  during  most  of  the  period  0  -  24  h. 
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6.  Forecasts  with  Assimilation  of  Low-level  Winds 

In  Exps.  4,  5,  and  6,  the  observed  low-level  radial  (ug)  and  tangential 
(vg)  winds  at  0  h  are  assimilated  by  DIR  into  the  24-h  forecase  during  a 
preforecase  integration  from  -12  to  0  h  (Table  2).  The  relaxation  coef¬ 
ficients  are  Xw,Xs,  and  in  Exps.  4,  5,  and  6,  respectively  (Figure  2). 

Figures  8  and  9  show  the  minimum  pressures  and  the  maximum  wind 
speeds  for  these  three  forecasts  with  low-level  winds  assimilated.  It  is 
apparent  that  the  DIR  with  the  weak  relaxation  coefficient  (Exp.  4)  does 
not  alter  the  prediction  appreciably  toward  the  observation  during  the  pre¬ 
forecast  integration.  The  following  24  h  forecast  shows  no  apparent  improve¬ 
ment  over  the  standard  forecasts. 

The  model  adjustments  are  considerable  when  strong  and  attenuating  relaxa¬ 
tion  coefficients  are  applied  in  Exps.  5  and  6.  The  maximum  wind  speed  in 
the  preforecast  integration  converges  to  the  observed  value  within  a  couple 
of  hours.  The  minimum  pressure  also  approaches  the  observed  value  at  0  h 
within  6  h,  in  agreement  with  the  theory  of  geostrophic  adjustment.  The 
maximum  wind  speed  in  Exp.  5  achieves  the  observed  value  due  to  constantly 
strong  relaxation.  However,  as  the  assimilation  terminates  at  0  h,  strong 
model  adjustments  occur  in  both  experiments.  The  model  rejection  occurs  in 
Exp.  5  where  the  maximum  wind  speed  decreases  approximately  5  m  s-1  in 
three  hours  and  stays  lower  than  that  of  the  standard  forecast  (Exp.  2).  The 
rejection  is  similar  in  Exp.  6,  so  there  is  still  no  improvement  in 
intensity  forecast.  ’ 

We  thus  conclude  that  low-level  wind  observations  are  not  beneficial 
to  intensity  forecasts  of  tropical  cyclones  if  assimilated  by  DIRT.  We 
can  also  conclude  that  the  "attenuating-'  relaxation  coefficient  X.  is  more 

a 

effective  in  assimilating  the  observed  data  (cf.  Fig.  9)  and 
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desirable  for  eliminating  model  adjustments. 

However,  an  examination  of  forecast  errors  in  Exps.  4,  5,  _  ___ 

and  6  is  warranted.  During  the  preforecast  integration,  e(V)  decreases 
with  time  as  relaxation  forces  the  low-level  wind  to  asymptotically  approach  the 
0  h  observation.  The  e(V)  value  from  Exp.  5  at  0  h  reaches  the  lowest  level 
of  all  (Figure  10).  The  model  adjustments  a  use  the  error  to  be  at  levels 
higher  than  those  of  the  standard  forecasts  (Exps.  2  and  3)  after  6  h  since 
the  assimilation  has  been  rejected  by  the  model.  The  e(T)  and  e(q)  are 
similar  to  e(V)  in  that  they  decrease  with  time  in  the  preforecast  integration 
for  strong  relaxation  and  they  subsequently  increase  to  levels  equivalent 
or  higher  than  those  of  the  standard  forecasts. 

The  rejection  of  assimilation  of  low-level  winds  in  the  preceding  experi¬ 
ments  can  be  attributed  to  the  insufficient  vertical  coupling  between  the  low- 
level  and  high-level  momentum  fields  during  the  12-h  period  of  preforecast 
integration.  A  longer  period  of  preforecast  integration  may  produce  enough 
vertical  coupling  through  model  dynamics  and  physics,'  but  is  not  very 
meaningful  in  practice.  It  is  then  logical  to  test  assimilation  of  additional 
wind  observations  at  higher  levels  since  they  can  be  made  available 
(Rodger  et  al^,  1979). 


7.  Forecasts  with  Assimilation  of  Low  and  Hi jher- level  Minds 

As  listed  in  Table  2,  observations  of  higher  level  winds  in  addition 
to  the  low-level  winds  are  assimilated  by  DIR  using  relaxation  coe- 
t  fficient  Xfl.  In  Exp.  7,  Vg  (check  Table  1  for  pressure  level)  at  0  h  is 

assimilated;  in  Exp.  8,  v^  and  v^;  and  in  Exp.  9,  U2  and  v£. 

The  assimilations  of  higher  level  winds  yield  significantly  different 

*  intensity  forecasts  from  the  standard  forecasts  as  evident  in  Figures  11 

and  12  showing  the  minimum  pressure  and  the  maximum  wind  speed,  respectively. 

In  addition,  the  three  experiments  forecast  very  different  minimum  pressures 

even  during  the  preforecast  integration  where  the  same  X  is  used.  Among  the 

3. 

three,  Exp.  8  yields  the  best  prediction,  with  maximum  difference  of  only  4.5 
mb  in  central  pressure  at  24  h.  Exp.  7  predicts  a  very  intense  storm,  with 
central  pressure  deepening  to  955  mb  at  12  h.  Exp.  9  predicts  a  weaker 
storm  than  the  observation,  however,  it  forecasts  better  than  the  standard 
forecasts  (Exps.  2  and  3) and  forecasts  where  only  low-level  winds  are 
assimilated  (Exps.  4,  5,  and  6). 

It  is  interesting  that  the  storm  intensities  in  Exps.  7  and  8  are 
drastically  different  when  the  only  difference  in  the  experiments  is  that 
the  observations  of  v^  at  0  h  are  available  for  assimilation  ir  Exp.  8.  As 
demonstrated  by  the  vertical  profile  of  v  at  r  =  30  km  in  Figure  13,  there 
is  a  strong  vertical  shear  in  the  0  h  tangential  wind  observation.  Tangen¬ 
tial  wind  speed  decreases  upward  associated  with  the  strong  warm  core  at 
level  4.  In  Exp.  7t  the  too  strong  storm  intensity  is  due  to  the  assimilation 
of  the  stronger  circulation  below  800  mb,  whereas  the  vertically  decreasing 
tangential  circulation  and  effects  of  warm  core  are  properly  assimilated  in 
Exp.  8.  This  suggests  that  when  observational  data  are  to  be  vertically 
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interpolated  in  diagnosis  or  analysis,  strong  vertical  shear  and  the 
related  baroclinic  effect  must  be  taken  into  account. 

The  forecast  errors  in  wind  speed,  temperature  and  water  vapor  are 
shown  in  Figures  14,  15,  and  16,  respectively.  The  values  of  e(V)  decrease 
in  time  during  the  preforecast  integration  as  in  Exps.  4,  5,  and  6.  Note 
the  e(V)  in  Exp.  9  is  the  smallest  because  initial  errors  are  largest  at 
the  assimilated  levels  (Figure  1).  The  errors  in  Exps.  8  and  9  remain 
smaller  than  those  of  the  standard  forecasts,  especially  in  Exp.  9,  where 
the  error  is  50 %  lower. 

The  error  in  temperature  field  of  Exp.  7  arises  early  in  the  prefore¬ 
cast  integration.  This  shows  that  the  effects  of  the  warm  core  on  tangential 
circulation  are  not  properly  assimilated  as  mentioned  earlier.  The  errors 
in  Exps.  8  and  9  are  generally  smaller  than  for  the  standard  forecasts 
throughout  the  24  h  forecast  period.  The  specific  humidity  errors  for  these 
three  forecast  experiments  are  higher  than  for  the  standard  forecasts  with 
Exp.  7  having  the  highest  error.  Between  -12  and  6  h,  e(q)  in  Exp.  8  is  very 
low  because  the' inflow  and  outflow,  which  nearly  determine  the  net  total 
water  vapor  convergence,  are  assimilated. 

The  higher  e(q)  in  Exps.  8  and  9, in  spite  of  the  better  intensity 
forecasts  and  lower  i(V)  and  e(T),  is  probably  due  to  the  different  physical 
parameterizations  in  the  forecast  model. 
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8.  Deterioration  of  Forecast  due  to  Satellite  Observation  Errors 

In  previous  sections,  we  have  examined  the  impact  of  assimilating  the 
error-free  satellite  winds.  The  assumption  of  perfect  satellite  observation 
is  for  clarity  in  comparison.  Satellite-sensed  winds  are,  of  course,  not 
error-free.  As  discussed  earlier,  the  mean  speed  error  in  satellite  winds 
varies  from  2.5  m  s~^  (Rodgers  et_  al_,  1979)  to  as  large  as  8  m  s-1  (Black, 
1979).  We  now  turn  our  attention  to  the  influence  of  the  satellite  observa¬ 
tion-errors  on  the  forecast. 

Because  low-  and  outflow-level  winds  are  most  likely  to  be  obtainable 
operationally ,  we  repeat  Exp.  9  with  artificially  introduced  observation 
errors  in  both  the  initial  condition  and  satellite-sensed  winds.  Randomized 
error  of  about  2.5  m  s"^  and  biased  errors  are  added  into  the  observation 
at  -12  h.  The  biased  errors  have  a  maximum  of  8  m  s~^  at  r  =  30  km, 
decreasing  with  radius  to  2ero  at  r  =  150  km.  These  errors,  after  the 
balanced,  static  initializations.,  are  equivalent  to  approximately  0.5°K 
errors  in  the  temperature  field  and  a  4  mb  error  in  central  pressure.  Four 
experiments,  Exps.  2E,  3E,  9E,  and  9E24are  carried  out  based  on  initial 
conditions  containing  such  errors.  The  satellite-sensed  winds  at  0  h  con¬ 
tain  random  errors  with  maximum  speed  error  of  2.5  m  s-^  in  Exp.  9E  and 
5  m  s  ^  in  Exp.  9E2.  Exps.  2E  and  3E  are  identical  to  Exps.  2  and  3  except 
for  the  introduced  errors  in  the  initial  condition. 

As  summarized  in  Table  3,  the  average  12- and  24-h  forecasts  in  Exps. 

2E,  3E,  and  9E  are  worse  than  their  respective  counterparts  in  error-free 
simulations.  For  example,  Exp.  9  has  an  averaged  forecast  error  of  11  mb 
in  minimum  pressure,  6.5  m  s~^  in  maximum  wind  speed,  and  2.5  m  s"^  in  e(V), 
whereas  Exp.  9E  has  an  averaged  forecast  error  of  16  mb,  11.7  m  s"\  and 
4.5  m  s  \  It  is  encouraging  that  Exp.  9E,  in  which  the  magnitude  of  the 
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Table  3.  12  and  24  h  Averaged  Forecast  Errors 


t 


t 


t 
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observation  errors  are  typical  for  operational  forecasts,  is  a  better 
forecast  than  the  standard  forecasts  of  Exps.  2E  and  3E.  It  forecasts 
better  than  the  standard  forecasts  by  4  mb  in  minimum  pressure,  approxi¬ 
mately  2  m  s'1  in  maximum  wind,  and  0.5  in  e(V).  It  is  also  interesting 
that  Exp.  9E2  performs  only  slightly  worse  than  Exp.  9E,  although  the 
error  level  is  twice  as  large  as  that  in  Exp.  9E. 

From  the  above  comparison,  we  conclude  that  the  errors  in  satellite 
winds  could  deteriorate  a  forecast,  and  that  the  assimilation  of 
satellite  low-  and  outflow-level  winds  can  improve  the  forecast  if  these 
errors  are  less  or  equal  to  those  contained  in  the  initial  wind  field. 
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9.  Summary  and  Discussion 

The  impact  of  accurately  measured  marine  surface  winds  of  sufficient 
spatial  coverage  and  resolution  in  the  24  h  intensity  forecast  of  tropical 
cyclones  has  been  studied  with  simulation  experiments.  The  model  physics 
in  the  forecast  model  were  altered  from  those  in  the  nature  model.  The 
■observations  are  assimilated  into  the  numerical  forecast  with  dynamical 
initialization  by  relaxation  during  the  preforecast  integrations  from 
error-free  mass  fields. 

The  results  indicate  no  improvement  in  forecast  accuracy  when 
low-level  winds  are  assimilated  according  to  the  abovementioned  pro¬ 
cedure.  We  note  that  a  strong  relaxation  coefficient  causes  rejection 
of  the  assimilation  within  a  few  hours  of  forecasting  and  that  a  weak 
relaxation  coefficient  is  ineffective. 

Significant  improvements  are  achieved  when  all  winds  below  600  mb 
are  assimilated.  This  conclusion  is  easily  understandable  because  the  major 
characteristics  of  the  tropical  cyclone  such  as  the  vortex  strength,  the 
warm  core,  and  the  vertical  shear  are  included  in  such  observations.  But 
simultaneous,  hfgh-resolution  observation  required  for  such  assimilation 

are  difficult  to  obtain.  It  is  encouraging  that  forecast  improvement 
can  also  be  achieved  when  low-and  outflow-level  winds  are  assimilated 
because  wind  fields  at  these  two  levels  are  most  likely  available  from 
satellite  observations. 

The  forecast  with  low-and  outflow-level  winds  assimilated  worsens  with 
increasing  observation  errors.  However,  even  if  the  root -mean-square 
error  in  the  satellite  observation  is  equivalent  to  that  in  the  initial 
wind  field,  assimilation  of  low-and  outflow-level  winds  still  improves  the 
forecast. 
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Caution  must  be  taken  in  interpreting  these  findings  for  operational 
I  applications,  as  is  the  case  for  all  simulation  studies,  because  the  extent 

to  which  they  approximate  reality  is  difficult  to  determine.  The  finding 
that  the  low-level  wind  observations  alone  cannot  improve  the  forecast  when 
f  assimilated  by  DIR  should  not  cast  doubt  on  the  usefulness  of  observing 

systems  which  measure  marine  surface  winds.  Since  an  axisymmetric  tropical 
cyclone  is  employed  in  this  study,  the  position  of  the  storm  is  assumed 
t  known.  Also,  the  mass  fields  are  assumed  to  be  error- free  in  the  static 

initializations  of  the  forecasts.  The  precise  center  location  and  perfect 
mass  field  are  not  commonly  available  for  operational  forecasts,  where 
*  meteorologists  have  to  be  content  with  uncertainties  of  the  storm  center 

position  and  with  the  "bogussed"  circulations.  The  marine  surface  winds 
are  invaluable  in  defining  the  low-level  circulations  and  in  locating  the 
storm  centers  which  otherwise  would  be  impossible  over  data-void  oceans. 

Since  our  results  with  assimilation  of  the  low  and  outflow-winds  are 
encouraging  and  since  these  will  be  the  focal  levels  in  satellite  observa¬ 
tions,  future  research  with  a  three-dimensional  tropical  cyclone  model  is 
warranted.  In  a  three-dimensional  study,  the  impact  of  the  satellite- 
sensed  winds  on  storm  track  forecasts  can  be  investigated.  The  effects  of 
time-lag  within  one  satellite  revolution  discussed  earlier  and  the  effects 
of  the  swath  width  can  also  be  studied.  Finally,  real  data  case  studies 
can  be  carried  out  with  a  three-dimensional  model. 
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ABSTRACT 


A  three-dimensional  numerical  model  with  a  domain  of  3000  x  3000  km  and 
horizontal  resolution  of  60  km  is  used  to  study  the  influence  of  sea- 
surface  temperature  (SST)  on  the  behavior  of  tropical  cyclones  translating 
with  mean  flows  in  the  Northern  Hemisphere. 

We  find  that  tropical  cyclones  tend  to  move  into  regions  of  warmer  SST 
when  a  gradient  of  SST  is  perpendicular  to  the  mean  ambient  flow  vector 
(MAFV).  The  model  results  also  indicated  that  a  region  of  warmer  SST  situated 
to  the  right  side  of  the  MAFV  is  more  favorable  for  storm  intensification  than 
to  the  left  side  due  to  the  asymmetries  in  air-sea  energy  exchanges  associated 
with  translating  tropical  cyclones.  The  model  tropical  cyclone  intensi¬ 

fies  and  has  greater  rightward  deflection  in  its  path  relative  to  the  MAFV 
when  translating  into  the  region  of  warmer  SST.  The  model  tropical  cyclone 
intensifies  when  its  center  travels  along  a  warm  strip,  it  weakens  along,  but 
does  not  move  away  from  a  cool  strip. 

The  results  suggest  that  the  SST  distribution  not  only  affects  the  intensity 

and  path  of  tropical  cyclones  fractionally,  but  also  affects  them  thermally. 

The  enhanced  evaporation  and  convergence  over  the  warm  SST  provide  a 
favorable  condition  for  the  growth  of  the  tropical  cyclone,  and  lead  to  a 
gradual  shift  of  the  storm  center  toward  the  warm  ocean. 


.1.  INTRODUCTION 


Unlike  their  counterparts  in  mid-  and  higher-latitude,  where  the  contrast 
of  air  masses  provides  enough  energy  for  circulation,  tripical  cyclones  de¬ 
pend  on  the  ocean  for  their  energy  source.  Sensible  and  latent  heat  trans¬ 
ported  through  the  air-sea  interface,  in  cooperation  wit)  the  frictionally- 
induced  inflow,  supplies  the  energy  for  maintaining  the  iroDical  cyclone 
circulation.  The  air-sea  interaction  is  an  important  physical  process  that 
affects  the  behavior  of  tropical  cyclones. 

Observational  studies  on  the  interaction  of  tropical  cyclones  and  the 
ocean  have  indicated  a  strong  dependence  of  tropical  cyclones  on  the  sea- 
surface  temperature  (SST).  The  development,  intensification,  and  movement  of 
tropical  cyclones  are  in  many  instances  linked  with  the  warm  SST  (Palmen, 

1948;  Fisher,  1958;  Miller,  1958;  Perlroth,  1967;  Brand,  1971;  Gray,  1979). 
Because  the  ocean  is  regarded  as  the  heat  reservoir,  attention  has  also  been 
focussed  on  the  importance  of  the  thermal  structure  of  the  upper  ocean.  It 
is  found  that  warm  and  deep  oceanic  mixed  layers  are  necessary  conditions  for 
the  development  and  intensification  of  tropical  cyclones  (Tisdale  and  Clapp, 
1963;  Perlroth,  1969).  Contrary  to  the  above  findings,  some  tropical  cyclones 
are  not  affected  by  SST  (Rcmage,  1972;  1974).  Because  many  other  factors 
can  influence  the  behavior  of  tropical  cyclones,  it  is  very  difficult  to 
isolate  the  effects  of  SST  from  observational  data. 

Earlier  numerical  studies  with  axi symmetrical  (therefore  stationary)  models 
showed  that  tropical  cyclones  are  very  sensitive  to  small  changes  of  SST  (Oo- 
yama,  1969,  Rosenthal,  1971;  Sundquist,  1972).  Typically,  a  50%  change  in 
model  maximum  wind  would  result  from  a  2°C  change  of  SST.  By  incorporating 
a  more  sophisticated  parameterization  of  the  atmospheric  boundary  layer  (ABL), 
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Anties  and  Chang  (1978)  showed  th  .  '  trie  tropical  cycle  no  is  not 

as  sensitive  as  earlier  studies  imiica'  d.  The  enhanced  (decreased)  evapora¬ 
tion  associated  with  warm  (cool)  SST  i  partially  compensatec  for  by  the 
increased  (decreased)  surface  friction,  therefore,  the  tropical  cyclone  has 
a  12h  delay  in  response  to  SST  changes.  When  their  tropical  cyclone  was 
coupled  with  the  ocean,  the  response  was  even  weaker  and  further  delayed 
(Chang  and  Anthes,  1979).  However,  the  long  term  behavior  of  their  model 
cyclone  still  agreed  with  previous  axisymmetric  model  studies,  i.e.,  warmer 
SST  results  in  intensification  and  cooler  SST  results  in  wealening. 

The  motion  of  a  tropical  cyclone  with  mean  flow  over  an  ocean  of  uni  form 
SST  has  been  studied  by  Kuo  (1969)  and  Jones  (1977a)  using  linear  models.  The 
air-sea  interactions  in  these  two  linear  models  were  purely  mechanical  with 
only  surface  friction.  These  studies  indicated  that  the  trajectory  of  the 
vortex  center  is  a  damped  trochoid.  The  deflection  angle  of  the  mean  path 
(the  angle  between  the  mean  current  and  the  mean  vortex  motion)  is  a  function 
of  the  surface  friction  and  the  environmental  friction.  For  tropical  cyclones, 
where  the  surface  friction  is  greater  than  the  environmental  friction,  the 
deflection  is  usually  to  the  right  of  the  mean  flow.  These  results  were  later 
substantiated  by  Jones’  (1977)  three-dimensional,  non-linear  model  with  uni¬ 
form  SST. 

The  behavior  of  tropical  cyclone  in  a  three-dimensional  model  with  mean 
flow  over  non-uniform  SST  has  not  been  studied.  With  variations  in  SST, 
both  the  intensity  and  movement  of  tropical  cyclones  may  be  more  complex. 

The  surface  friction  will  change  because  of  variationsin  surface  stability 
associated  with  the  non-uniform  SST,  the  sensible  and  latent  heat  exchanges 
may  be  in  certain  ways  dictated  by  the  SST  distributions.  These  variations 
may  significantly  affect  the  behavior  of  tropical  cyclones. 
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The  purpose  of  this  study  is  to  investigate  the  influence  of  var¬ 
ious  non-uni  form  SST  distributions  on  the  intensities  and  movements  of  tror  - 
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ical  cyclones  that  are  imbedded  in  a  mean  flow.  As  described  in  the  next 
section,  the  model  domain  is  large  enough  to  allow  free  movement  of  the  trcD- 
ical  cyclone  for  about  3-4  days.  The  parameterization  of  the  ABL  has  been 
improved  so  that  the  air-sea  interaction  is  adequately  treated.  In  Sectior 
III,  the  life  history  of  a  weak  tropical  cyclone  translating  over  ocean  with 
uniform  SST  will  be  discussed  and  used  as  the  control  experiment  for  later 
comparison.  In  Section  IV,  we  first  examine  the  responses  of  tropical  cy¬ 
clones  when  translating  into  regions  of  warmer  and  cooler  SST.  In  Section 
V,  the  response  of  the  translating  tropical  cyclone  with  SST  gradients  perpen¬ 
dicular  to  the  mean  flow  will  be  discussed.  Finally,  the  repsonse  of  the 
translating  tropical  cyclone  encountering  strips  of  warmer  or  cooler  SST  will 
be  discussed  in  Section  VI. 

II.  THE  NUMERICAL  MODEL 

The  tropical  cyclone  model  used  is  identical  to  the  one  described  in 

(1) 

Madala  and  Chang  (1979).  The  governing  equations  include  the  primitive  equa¬ 
tions  of  conservation  for  horizontal  momentum,  mass,  enthalpy,  and  water  vapor. 
The  system  of  equations  is  hydrostatic.  A  normalized  pressure  (o)  is  the  verti¬ 
cal  coordinate  (Phillips,  1957).  The  physics  includes  the  subgrid  scale  hori¬ 
zontal  mixing,  the  subgrid  scale  cumulus  convection,  the  grid  scale  precipita¬ 
tion,  and  the  subgrid  scale  vertical  mixing  due  to  surface  friction. 

The  subgrid  scale  horizontal  mixing  is  parameterized  by  a  kinematic  eddy 
coefficient.  This  coefficient  consists  of  a  constant  part  and  a  part  linearly 
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dependent  on  wind  speed.  This  form  of  the  eddy  coefficient  yields  suitable 
mixing  in  the  initial  as  well  as  the  mature  stages  (Anthes  et  al^  1971). 

The  subgrid  scale  cumulus  convection  is  parameterized  following  Kuo's 
(1974)  and  Anthes's  (1977)  methods.  Conditional  instability  and  the  boundary 
layer  convergence  of  water  vapor  are  prerequisites  for  the  cumulus  convection. 
The  partitioning  of  heating  and  moistening  depends  on  the  vertically  averaged 
relative  humidity  of  the  air  column.  The  vertical  distribution  of  heating  is 
determined  by  the  amount  of  conditional  instability  and  the  prescribed  weight¬ 
ing  functions  (see  Anthes  and  Chang,  1979).  The  use  of  the  weighting  functions 
insures  a  proper  vertical  heating  distribution  for  the  growth  of  tropical  cy¬ 
clone  scale  disturbances.  The  grid  scale  precipitation  occurs  when  the  air 
reaches  saturation  in  grid  scale  lifting. 

The  boundary  layer  effects  are  parameterized  based  on  a  generalized  sim¬ 
ilarity  theory  in  which  the  logarithmic-linear  profiles  of  the  momentum, temp¬ 
erature,  and  water  vapor  are  "match"  into  the  mixed  layer  (Chang  and  Madala, 

•  (2) 

19  80).  Anthes  and  Chang  (1978)  have  shown  a  higher  resolution  BL  is  superior 
to  a  single  layer  bulk  ABL  in  an  axisymmetric  tropical  cyclone  model.  To 
incorporate  into  a  3D  tropical  cyclone  model  a  higher  resolution  ABL  is  econo¬ 
mically  prohibitive.  In  addition,  the  results  of  Anthes  and  Chang  (1979)  indi¬ 
cate  that,  except  for  the  period  of  rapidly  changing  stability,  the  ABL  in 
tropical  cyclones  with  strong  winds  and  small  Richardson  numbers  is  dominated 
by  mechanical  mixing  and  is  nearly  invariate  in  height.  Therefore,  matching 
technique  with  fixed  ABL  height  seems  to  be  appropriate.  Universal  functions 
A,  B,  C,  and  D  are  formulated  following  Yamada  (1976).  The  surface  roughness 
height  over  open  water  are  calculated  by  Charnock's  equation  (Del sol  et  ajh 


The  model  <:mosphere  from  P  =  P:  (surface)  to  P  =  0  is  divided  into 
seven  sigma  layers  (Figure  1).  Note  lhat  the  lowest  layer  has  a  thickness 
of  about  700  m,  which  is  the  typical  /  BL  depth  in  tropical  cyclones  (Moss 
and  Merceret,  1976;  Anthes  and  Chang,  1978;  Chang  and  Anthes,  1979).  A 
better  vertical  resolution  between  o  -  0.1  and  o  =  0.3  is  used  to  resolve  the 
tropopause  and  outflow  structure.  .  All  prognostic  variables  such  a 
and  q,  are  defined  at  the  center  of  e;  ch  layer,  all  diagnostic  variables 
such  as  o  and  u>  are  defined  at  the  bot ndary  of  each  layer.  The  momentum  points 
and  mass  points  are  fully  staggered  ir  horizontal  direction  following  the 
scheme  C  in  Arakawa  and  Lamb  (1977).  The  uniform  horizontal  resolution  is  60 
km  throughout  the  model  domain  of  300C  x  3000  km.  We  take  note  that  a  hori¬ 
zontal  resolution  of  60  km  is  marginal  in  resolving  the  detail  structure  near 
the  storm  center.  However,  a  resolution  of  60  km  is  a  reasonable  compromise 
for  the  time  and  space  scales  of  this  study.  The  spatial  differencing  is  second 
order  to  conserve  mass,  momentum,  and  enthalpy.  The  newly  developed  split- 
explicit  method  is  applied  for  temporal  integration.  The  boundary  conditions 
are  Neuman  for  all  dependent  variables.  A  sponge  boundary  condition  is 
applied  to  the  momentum  to  control  noise  at  lateral  boundaries. 

The  tropical  atmosphere  of  the  hurricane  season  (Sheets,  1971)  is  used 
for  the  initial  mean  thermodynamic  state  of  the  model.  The  thermal  structure 
features  a  conditionally  unstable  lapse  rate  from  surface  to  about  350  mb  and 
high  relative  humidity  (RH)  up  to  500  mb. 

The  initial  flow  field  includes  a  uniform  easterly  current  of  about 

4.5  m  s'1  and  a  non-divergent,  idealized  vortex  with  maximum  wind  of  16  m  s  . 

The  initial  surface  pressure  and  temperature  field  are  in  gradient  balance 

-5  -1 

with  the  flow  field.  The  Coriolis  parameter  is  5  x  10  s  .  One  effect 


of  latitudinal  variation  of  Coriolis  para:  •.<.  -_er  is  to  produce  a  i.orthward 
movement  of  the  vortex  (Madala  and  Piacset. ,  1975;  Anthes  and  Ho/.e,  1975). 
For  the  purpose  of  isolating  the  movement  of  tropical  cyclones  ;aused  by 
SST ,  the  Coriolis  parameter  is  set  constant  in  this  study. 
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III.  THE  CONTROL  EXPERIMENT 


A  tropical  cyclone  with  rather  weak  intensity  develops  from  the  initial 
vortex  imbedded  in  the  uniform  easterly  mean  flow  of  approximately  4.5  m  s~^ 
over  an  ocean  of  28.5  C.  The  model  tropical  cyclone  reaches  a  quasi -steady 
state  after  48  h,  attaining  a  minimum  pressure  of  996  mb~^  and  maximum  ABL 
wind  speed  in  the  lowest  layer  of  the  model  of  32  m  s"^  (this  experiment 
is  referred  to  as  the  control  experiment  (Exp.  A,  see  Table  1)). 

Figure  2  shows  the  ABL  flow  pattern  and  the  surface  pressure  field 
at  84  h  when  the  model  integration  was  terminated.  The  strong  circulation 
is  confined  to  the  region  near  the  storm  center.  Maximum  winds  are  located 
to  the  right  of  the  storm  track  where  the  cyclonic  circulation  and  the  mean 
current  is  in  the  same  direction.  The  winds  at  the  center  are  very  weak. 

The  flow  far  from  the  storm  circulation  is  nearly  undisturbed  by  the  storm. 
These  characteristics  of  the  low-level  flow  are  also  implied  by  the  surface 
pressure  field.  The  circular  isobars  and  the  stronger  pressure  gradient  are 
concentrated  to  the  region  near  the  storm  center. 

Contrary  to  the  concentrated  cyclonic  circulation  at  low  level,  the 
asymmetric  flow  at  the  outflow  level  (a  =  0.1  -  0.2)  is  predominantly  anti- 
cyclonic  with  a  small  region  of  cyclonic  circulation  near  the  center.  The 
speed  of  the  anticyclonic  circulation  increases  away  from  the  center  in 
conservation  of  the  angular  momentum.  It  reaches  a  maximum  of  more  than 
20  m  s'1  to  the  southeast*  of  the  storm  center. 


Because  the  Coriolis  parameter  is  constant,  the  orientation  of  the  domain 
has  no  real  geophysical  meaning. 
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As  shown  by  the  paths  of  the  sto^m  center  (defined  by  the  surface 
pressure  field),  the  storm  moves  downstream  with  the  mean  current  and  with 
a  deflection  to  the  right  of  the  MAFV.  At  84  h,  the  net  northward  displace¬ 
ment  of  the  storm  center  is  about  100  km.  The  5°  deflection  angle  of  our 
model  is  very  typical  for  model  tropical  cyclones  (Jones,  1977a).  However, 
the  small  scale  looping  motion  in  the  trochoid  (Jones,  1977b)  cannot  be 
resolved  by  the  60  km  grid  resolution  of  the  present  model. 

Figure  3  shows  the  stress  field  at  84  h  of  Exp.  A.  Comparing  with 
the  surface  wind  field,  the  region  of  the  strongest  stress  is  located  to 
the  right  of  the  storm  track.  The  crescent  shape  of  Lne  maximum  stress 
region  is  responsible  for  the  often  observed  crescent-shaped  maximum  cooling 
in  the  ocean  (Black,  1972).  To  the  south  of  the  strong  stress,  there  is  a 
region  where  the  stress  is  very  small,  corresponding  to  the  weak  surface 
wind  at  the  storm  center. 

The  evaporation  rate  at  84  h  of  Exp.  A  shows  that  the  strongest 
evaporation  occurs  outside  and  to  the  north  of  the  region  of  the  maximum 
wind  (Fig.  4)  where  the  ABL  winds  are  considerable  and  the  humidity  is 
undersaturated.  The  maximum  evaporation  rate  is  approximately  2.5  cm  day”\ 
The  sensible  heating  is  relatively  insigni ficant,  in  agreement  with  earlier 
findings  of  Rosenthal  (1971)  and  Anthes  and  Chang  (1978). 

A  series  of  experiments  with  various  non-uniform  SST  distributions 
were  performed.  The  schematical  patterns  of  SST  tested  are  shown  in 
Figures  5  and  listed  in  Table  1. 

IV.  TROPICAL  CYCLONES  TRANSLATING  INTO  BROAD  REGIONS  OF  WARM  AND  COOL  SST 

As  mentioned  before,  the  intensity  of  axisymmetric,  stationary  model 
tropical  cyclones  is  sensitive  to  sudden  changes  of  SST.  The  situation  that 
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Table  1.  SUMMARY  OF  EXPERIMENTS  AND  RESULTS 


we  consider  here  is  a  tropical  cyclone  that  translates  into  regions  of 
warmer  or  cooler  SST.  The  response  of  the  tropical  cyclone  in  this  more 
realistic  situation  is  investigated  with  Exp.  W. ,  where  the  SST  is  2  C 
warmer  at  30.5  C  in  the  western  half  of  the  domain,  and  Exp.  C,  where  the 
SST  is  2  C  cooler  at  26.5  C  in  the  western  half  of  the  domain.  SST  in  the 
eastern  half  of  the  domain  stays  at  28.5  C. 

Figure  6  shows  the  time  series  of  the  minimum  pressures  of  Exps.  W  and 
C  and  the  control  (Exp.  A).  In  both  Exps.  W  and  C,  the  minimum  pressures 
deviate  very  little  from  the  control  experiment  until  the  center  of  the 
storm  moves  into  the  region  of  warmer  or  cooler  SST.  The  delay  is  also 
evident  in  the  maximum  surface  wind  which  does  not  deviate  significantly  from 
the  control  until  after  36  h  (Figure  7).  This  delayed  response,  similar  to 
the  behavior  of  the  axi symmetric  model  of  Anthes  and  Chang  (1978),  is  due  to 
the  increased  (decreased)  dissipation  of  kinetic  energy  which  partially 
compensates  for  the  enhanced  (weakened)  evaporation  over  the  warm  (cool)  SST. 

The  SST  patterns  tested  here  also  affect  the  path  of  the  tropical  cyclone. 
At  84  h,  the  storm  centers  in  Exps.  W  and  C  are  displaced  to  about  60-100  km 
to  the  right  and  left,  respectively,  relative  to  the  control  experiment  (Fig¬ 
ure  8).  The  fact  that  the  increased  surface  friction  produces  a  stronger  de¬ 
flection  to  the  right  of  the  MAFV  agrees  qualitatively  with  the  linear  models 
of  Kuo’s  (1969)  and  Jones'  (1977a),  But  in  contrast  to  the  linear  theory,  the 
variation  in  the  surface  friction  may  be  due  to  both  storm  intensity  and  the 
surface  stability.  As  it  becomes  clearer  in  the  experiments  discussed  in  the 
next  section,  the  asymmetries  associated  with  a  translating  tropical  cycline 
and  air-sea  thermal  exchange  can  also  contribute  to  the  deflection.  Most  of 
the  air-sea  exchanges  of  sensible  and  latent  heats  occur  in  the  right  half  of 
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the  tropical  cyclone  relative  to  the  MAFV,  and  have  important  thermal 
effects  on  the  paths  of  tropical  cyclones  other  than  the  mechanical  effects 
considered  previously.  As  shown  in  Figure  9,  stronger  evaporation  mainly 
occurs  over  the  region  of  warmer  SST  at  27  h  of  Exp.  W.  The  maximum  evapora¬ 
tion  rate  increases  from  about  2  cm  day-^  at  24  h  to  more  than  6  cm  day”^ 
at  36  h  in  Exp.  W  when  the  center  of  the  storm  moves  into  the  region  of 
warmer  SST.  Such  a  strong  and  asymmetric  evaporation  distribution  causes  a 
stronger  convection  and  heating  to  the  right  of  the  MAFV  and  shifts  the 
storm  center  toward  the  north. 

V.  TROPICAL  CYCLONE  TRANSLATING  PERPENDICULAR  TO  SST  GRADIENT 

To  further  differentiate  ocean's  thermal  effect  on  the  path  of  tropi¬ 
cal  cyclones,  we  now  consider  the  case  that  the  mean  flow  is  perpendicular 
to  SST  gradient.  In  Exp.  WR,  the  SST  is  30.5  C  to  the  right  and  26.5  C  to 
the  left  of  the  MAFV,  and  in  between  the  30.5  C  and  26.5  C  waters,  there 
is  a  600-km  wide  region  where  the  SST  increases  continuously  toward  the 
right  having  a  gradient  of  4  C/600  km.  In  Exp.  WL,  the  SST  pattern  is 
reversed  so  that  the  warm  and  cool  ocean  waters  are  located  to  the  left  and 
right  of  the  MAFV,  respectively,  and  the  SST  gradient  is  toward  the  left 
i  n  the  central^  j^aianu-  — - 

These  SST  patterns  have  very  noticeable  influence  on  both  the  inten¬ 
sity  and  the  path  of  the  model  tropical  cyclone.  As  illustrated  in  Figs. 

10  and  11,  the  storm  intensity  of  Exp.  WR  increases  greatly  to  a  central 
pressure  of  982  mb  and  maximum  wind  of  50  m  s“^  at  84  h;  whereas  the 
storm  intensity  in  Exp.  WL  remains  essentially  unchanged.  The  storm  centers 


in  both  Exps.  WR  and  WL  are  displaced  perpendicular  to  the  MAFV  from  the 
control  experiment  by  more  than  200  km  toward  the  regions  of  warmer  SST 
(Fig.  12).  At  84  h,  the  center  of  circulation  in  Exp.  WR  has  almost  entered 
the  region  of  30. 5C  water.  The  apparent  tendency  of  the  model  tropical 
cyclone  to  move  toward  the  warm  SST  is  caused  by  the  enhanced  latent  heating 
associated  with  the  SST  distributions. 

As  illustrated  by  Figures  13  and  14,  the  evaporation  in  Exp.  WR  occurs 
mainly  over  the  warm  water  to  the  right  of  the  MAFV  and  to  the  rear  of  the 
storm  where  the  air  flows  into  warmer  water.  The  strong  evaporation  in  Exp. 
WL  is  obviously  biased  toward  the  southwest  of  the  storm  center.  Both 
the  rate  of  evaporation  and  the  area  of  strong  evaporation  in  Exp.  WR  are 
much  larger  than  in  Exp.  WL.  Due  to  the  asymmetry  associated  with  the 
translating  tropical  cyclone,  the  surface  wind  and  the  surface  stress  are 
stronger  to  the  right  of  the  MAFV  (see  Figure  2  and  3).  In  Exp.  WR  the 
region  of  stronger  stress  coincides  with  the  region  of  warmer  SST,  whereas 
in  Exp.  WL,  the  region  of  stronger  stress  does  not  coincide  with  the  region 
of  warmer  SST.  Consequently,  the  maximum  evaporation  rate  increases  to  6~7 
cm  d”l  in  Exp.  WR,  it  remains  at  about  2  cm  d~^  in  Exp.  WL. 

The  overall  stronger  evaporation  in  Exp.  WR  causes  the  storm  to  inten¬ 
sify.  Increased  evaporation  and  enhanced  surface  friction  also  cause  a  local 
increase  in  the  latent  energy  convergence  over  the  region  of  warmer  SST. 

Thus  provides  a  favorable  condition  for  the  growth  of  the  tropical  cyclone. 
The  stronger  convergence  and  hence  the  enhanced  cumulus  heating  over  the 
wanner  SST  cause  the  center  of  circulation  to  shift  toward  the  warm  water. 

The  shift  is  then  followed  by  the  pressure  field  through  adjustment  process. 
This  can  be  illustrated  by  the  isobar  patterns  and  flow  fields  at  48  h 
(Figures  13  and  14).  The  flow  fields  shown  are  at  the  level  of  o  =  0.66, 
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which  would  nearly  be  in  gradient  balance  with  the  pressure  field  under 
steady  state  and  stationary  conditions.  It  is  apparent  in  Figures  13  and  14 
that  the  circulation  centers  are  ahead  of  the  low  pressure  center  in  their 
movements  toward  the  warmer  SST.  The  isobars  are  not  concentric  with  weaker 
pressure  gradients  located  toward  the  warm  SST,  indicating  the  shift  of  the 

$ 

mass  field  toward  the  warmer  ocean. 

Note  also  that  the  rightward  deflection  from  the  control  path  in 
®  Exp.  WR  is  larger  than  that  in  Exp.  W,  and  the  leftward  deflection  from 

the  control  path  in  Exp.  WL  is  larger  than  that  in  Exp.  C,  where  the 
storm  intensity  is  weaker  than  Exp.  WL.  This  behavior  cannot  be  explained 
by  the  linear  theories  where  the  hurricane-ocean  interaction  is  merely 
mechanical.  The  SST  distribution  causes  a  redistribution  of  the  latent 
heating,  and  shifts  the  circulation  center  toward  the  region  where  the  growth 
of  the  tropical  cyclone  is  favorite. 

3 

Several  investigators  have  used  the  total  thermal  potential  as  a 
measure  of  the  energy  available  for  the  development  of  incipient  tropical 
cyclones  (Perlroth,  1969;  Leipper  and  Volgenan,  1972).  As 
illustrated  by  Exps.  WR  and  WL,  the  storm  intensities  and  paths  are  very 
different  in  spite  of  the  fact  that  the  total  oceanic  thermal  potentials 
in  two  cases  are  identical.  This  confirms  the  finding  in  Chang  (1979) 
that  the  air-sea  energy  exchange  in  the  tropical  cyclone  and  ocean  system 
depends  very  much  on  the  details  of  the  relative  location,  size,  and 
magnitude  in  SST  variations. 
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VI.  TROPICAL  CYCLONE  TRANSLATING  INTO  STRIPS  OF  WARM  OR  COOL  SST 


We  now  examine  the  response  of  the  tropical  cyclone  when  it  passes 
over  strips  of  warmer  or  cooler  SSTs.  As  depicted  in  Figure  5,  the  SST 
patterns  include  (1)  a  300-km-wide  2C  warm  strip  parallel  to  the  MAFV 
(Exp.  WS),  (2)  a  300-km-wide  2C  cold  tongue  parallel  to  the  MAFV  (Exp. 

CS),  and  (3)  a  300-km-wide  2C  warm  strip  at  45°  angle  to  the  right  of 
the  MAFV. 

The  response  of  the  nodel  tropical  cyclone  to  the  warm  or  cold  strip 
is  not  noticeable  until  after  36  h  when  the  center  of  the  storm  moves  in¬ 
to  the  warm  or  cold  tongue.  At  84  h,  the  minimum  pressures  of  Exps.  WS 
and  CS  are  respectively  980  and  1000  mb,  and  the  maximum  surface  winds  are 
respectively  50  and  25  m  s_1  (Figures  15  and  16). 

Contrary  to  the  finding  that  tropical  cyclones  tend  to  steer  away 
from  a  pool  of  cold  ocean  water  (e.g.,  Brand,  1971),  the  paths  of  the  storm 
centers  in  Exps.  WS  and  CS  are  essentially  unchanged  as  compared  to  the 
control  Exp.  A  despite  the  difference  in  storm  intensity  (Figure  13).  It  is 
probably  due  to  the  fact  that  the  changes  in  surface  friction  are  restricted 
to  the  relatively  small  regions  of  warm  or  cold  strips.  The  gradual  shift 
of  the  vortex  toward  the  region  of  warm  ocean  evident  in  Exps.  WR  and  WL  does 
not  occur  in  Exp.  CS  because  the  warmer  water  is  not  immediately  adjacent  to 
the  center  of  the  storm. 

The  model  tropical  cyclone  has  very  little  response  to  the  warm  strip 
which  lies  at  45  angle  to  the  MAFV  (Exp.  WS45).  Previous  experiments  with 
axi symmetric  hurricane  model  (Anthes  and  Chang),  1978;  Chang  and  Anthes, 

1979)  showed  that  a  tropical  cyclone  does  not  respond  to  warm  SST  if  the 
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tropical  cyclone  overlays  the  region  of  warmer  SST  for  less  than  12  h. 

*  The  storm  in  Exp.  WS45  travels  through  the  region  of  warmer  SST  in  less 
than  15  h,  which  is  apparently  not  long  enough  for  the  tropical  cyclone 
to  have  significant  response  before  it  leaves  the  warm  strip. 

! 

The  variation  in  SST  affects  the  air-sea  exchange  of  sensible  and 
latent  heat.  As  shown  by  Figure  14,  the  strong  evaporation  at  30  h  of 
I  Exp.  WS  occurs  mainly  in  the  region  of  the  warm  SST.  The  strongest  evap¬ 

oration  occurs  near  the  leading  (northern)  edge  of  the  warm  tongue.  The 
evaporation  in  Exp.  CS  is  much  weaker  over  the  cold  tongue  except  for  a 
>  slight  increase  where  the  air  flow  re-enters  the  region  of  28. 5C  water. 

Changes  in  the  evaporation  near  the  storm  center,  over  a  period  of  time, 
results  in  changes  in  the  intensities  of  the  tropical  cyclone.  There  are 
many  other  fine  details  in  the  momentum  field  of  the  tropical  cyclone  in¬ 
duced  by  the  warm  and  cold  strips.  For  example,  before  the  center  of  the 
tropical  cyclone  and  the  major  portion  of  the  strong  cyclonic  circulation 
'  enter  the  region  of  change  SST,  the  low-level  flow  reacts  to  the  warm  strip 

in  a  fashion  similar  to  that  of  a  heat  island.  Upward  motion  is  induced  in 
the  ABL  due  to  temperature  gradient  within  the  edge  of  the  warm  ocean. 

*  These  changes  in  the  momentum  field  are  persistent,  but  localized,  and  do 
not  seem  to  affect  the  overall  behavior  of  the  tropical  cyclone. 

VII.  SUMMARY 

The  influence  of  various  SST  distributions  on  the  intensities  and 
paths  of  translating  tropical  cyclones  has  been  investigated  using  a  three 
dimensional  model  of  tropical  cyclones  with  improved  parameterization  of  the 
ABL  and  Kuo's  (1974)  cumulus  parameterization. 
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Various  distribution  of  the  SST  are  tested  and  found  to  have  consider- 
able  influence  on  the  intensity  and  path  of  the  model  tropical  cyclone  (see 
Table  1).  The  variations  in  the  SST  not  only  changes  the  total  surface 
friction,  which  controls  the  deflection  of  the  vortex  from  the  mean  ambient 
flow,  but  also  alters  the  sensible  and  latent  heat  exchanges.  The  magnitude 
and  the  spatial  distribution  of  the  heat  exchanges  affect  the  intensity  as 
well  as  the  movement  of  the  trooical  cyclone.  Major  results  are  summarized 
as  follows: 

1.  There  is  a  deflection  of  approximately  five  degrees  to  the  right 
of  MAFV  as  a  slowly  varying  tropical  cyclone  translates  on  an  f-plane  with 

i 

mean  flow  over  an  ocean  of  uniform  SST. 

2.  The  intensity  and  the  angle  of  deflection  increases  with  increased 
SST  due  to  enhanced  evaporation  and  friction. 

3.  When  translating  downwind  the  mean  flow  over  the  ocean  with  SST 
gradients  perpendicular  to  the  MAFV,  tropical  cyclones  tend  to  move  into  the 
region  of  warmer  SST.  The  movement  toward  warm  SST  is  gradual  and  con- 
tinouous. 


4.  Tropical  cyclones  are  more  likely  to  intensify  when  warmer  ocean 
situates  to  the  right  (left)  of  the  storm  track  than  to  the  left  (right) 
in  the  Northern  (Southern)  Hemisphere  because  of  the  asymmetries  associated 
with  translating  tropical  cyclones. 
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5.  Narrow  (as  compared  to  the  storm  size)  regions  of  warmer  or 
pooler  SST  have  little  influence  on  the  paths  of  tropical  cyclones.  How¬ 
ever,  if  the  central  portion  of  the  tropical  cyclone  overlays  the  narrow 
region  of  warm  (cooler)  SST  over  a  period  of  time  (>  12  h),  the  tropical 
cyclone  intensifies  (weakens). 

From  these  results,  we  conclude  that  the  SST  distribution  affects  the 
behavior  of  a  translating  tropical  cyclone,  both  frictionally  and  thermally. 
In  agreement  with  the  adiabatic  linear  theories,  the  translating  tropical 
cyclone  appears  to  deflect  more  to  the  right  of  the  mean  flow  when  friction 
is  increased  (either  due  to  increased  surface  instability  or  stronger  in¬ 
tensity).  The  SST  distribution  also  redistributes  the  available  latent 
energy.  This  causes  a  shift  in  cumulus  heating,  provides  a  favorable  condi¬ 
tion  for  the  growth  of  the  tropical  cyclone,  and  ultimately  shifts  the 
vortex  toward  the  warm  ocean. 

In  light  of  these  results,  an  increased  effort  to  obtain  and  utilize 
accurate  SST  data  in  operational  forecast  of  tropical  cyclones  deserves 
more  serious  consideration. 
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1.  Madala,  R.  V.  ,  and  S.  W.  Chang,  1979:  A  vectorized  three-dimensional 

operational  tropical  cyclone  model.  Preprint  Scientific  Computer 
Information  Exchange  Conference,  August  1979,  Livermore. 

2.  Chang,  S.  W.  and  R.  V.  Madala,  1980:  Planetary  boundary  layer 

parameterization  for  tropical  cyclones  based  on  generalized  similarity 
theory.  Naval  Research  Laboratory  Memorandum  Report  4235.  Available 
at  Code  4780,  NHL,  Washington,  DC  20375. 

3.  The  areal  sum  of  the  vertically  integrated  heat  content  above  some 

critical  temperature. 
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ABSTRACT 


A  planetary  boundary  layer  (PBL)  parameterization  based  on  the 
generalized  similarity  theory  (GST)  was  tested  in  tropical  cyclone 
models.  This  parameterization,  with  only  one  layer,  is  desired  in 
modeling  tropical  cyclones  for  computational  speed.  The  momentum, 
sensible  heat  and  moisture  fluxes  are  mutually  dependent  in  this  para¬ 
meterization  through  non-dimensional  gradient  equations.  The  internal 
structure  of  the  PBL  is  determined  implicitly  through  universal 
functions. 

In  comparison  with  a  complex,  one-dimensional,  multi-layer  PBL 
model,  the  GST  parameterization  yields  accurate  moisture  fluxes,  but 
slightly  overestimates  the  momentum  flux  and  underestimates  the 
sensible  heat  flux.  The  GST  parameterization  produces  very  realistic 
dynamics,  energetics  and  thermal  structure  in  an  axi symmetric  tropical 
cyclone  model.  This  GST  parameterization,  although  unable  to  treat 
the  diffusion  across  the  PBL  inversion,  is  judged  superior  to  drag 
coefficient  parameterization  and  is  a  good  alternative  to  the  more 
expensive,  multi-layer  parameterization. 


I. 


INTRODUCTION 


The  planetary  boundary  layer  (PBL)  plays  a  critical  role  in 
the  evolution  of  tropical  cyclones  because  the  air-sea  energy  and 
momentum  exchanges  occur  through  the  PBL.  Parameterization  of  the 
PBL  is  thus  a  very  important  aspect  of  tropical  cyclone  modeling. 
Complex  parameterization  is  not  always  computationally  feasible  for 
three-dimensional  or  operational  models.  Computationally  efficient 
yet  accurate  PBL  parameterization  is  imperative  for  these  models. 

The  purpose  of  this  paper  is  to  test  such  a  parameterization  based  on 
the  generalized  similarity  theory  (GST)  in  the  parameterization  of  the 
PBL  in  tropical  cyclone  models. 

The  most  widely  used  PBL  model  is  the  single-layer  parameteri¬ 
zation  where  the  PBL  is  represented  by  one  model  layer  and  the  surface 
fluxes  are  calculated  using  bulk  aerodynamic  formulas.  This  method 
is  simple  and  computationally  economical  but  has  a  major  drawback  in 
that  the  internal  structure  of  the  PBL  cannot  be  resolved.  As  a  con¬ 
sequence,  the  single  layer  PBL  model  using  bulk  aerodynamic  formulas 
can  only  crudely  estimate  various  surface  fluxes. 

An  alternative  approach  where  the  PBL  is  modeled  by  several 
layers  has  been  adopted  in  recent  years  by  Pielke  (1974),  Kurihara 
and  Tuleya  (1974),  Anthes  and  Chang  (1978)  and  Anthes  and  Warner  (1978). 
By  explicitly  resolving  the  PBL  structure,  these  models  often  yield 
more  realistic  estimates  of  the  various  fluxes,  but  the  computational 
cost  is  high. 
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Recent  observational  studies  have  facilitated  the  testing  of 
various  PBL  similarity  theories.  The  rather  idealized  assumptions 
characteristic  of  early  theories  have  been  relaxed  and  improved,  and 


more  realistic  and  general  theories  have  been  derived  (Arya  1977; 

Yamada,  1976).  These  developments  have  led  to  a  PBL  formulation 
that  appears  suitable  to  model  various  PBL  conditions. 

In  the  next  section,  a  brief  review  of  the  generalized  similarity 
theory  will  be  given.  In  Section  3,  various  surface  fluxes  in  typical 
tropical  cyclone  conditions  obtained  from  the  GST  will  be  compared  with 
those  computed  by  a  time  dependent  one  dimensional  PBL  model  (Busch, 
et  al . ,  1976;  hereafter  referred  as  BCA).  The  generalized  similarity 
theory  is  then  applied  to  an  axi symmetric  tropical  cyclone  model. 

The  results  will  be  compared  in  Section  4  with  those  obtained  using 
constant  drag  coefficients. 


II. 


GENERALIZED  SIMILARITY  THEORY 

Similarity  theory  states  that  the  profile  of  any  properly  scaled 


variable  in  the  PBL  such  as  wind,  temperature,  or  water  vapor,  can  be 
described  by  a  universal  function.  Although  there  has  been  much  debate 
concerning  the  proper  scales  for  the  various  types  of  PBL's  (Arya,  1977) 
we  will  use  a  similarity  theory  that  appears  most  applicable  and  general 
for  tropical  cyclones  (Yamada,  1976). 

The  mathematical  representation  of  the  PBL  flow  can  be  obtained 
by  assuming  that  there  exists  a  layer  where  both  the  equations  de¬ 
scribing  the  surface  layer  and  the  equations  describing  the  interior 
flow  in  the  mixed  layer  are  both  valid.  In  the  surface  layer,  the 
PBL  variables  are  assumed  to  follow  a  logarithmic-linear  relationship 
with  height: 


|ln  (§-)  -  Tm< 


r> 


ill  .if, 

U*  k  L1 

r-  ■  [r  -  'H  <f»] 
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where  t  is  the  unit  vector  parallel  to  the  surface  wind,  Pr  is  the 
turbulent  Prandtl  number,  4^,  4'h,  and  V  are  stability  functions, 
zQ  is  the  roughness  length,  and  0Q  and  qQ  are  potential  temperature 
and  specific  humidity  at  zQ. 
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For  the  interior  flow  in  the  mixed  layer  we  can  write  the 


equations  in  the  form  of  the  resistance  law: 


V-V 
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where  V^,  6m  and  qm  are  the  mean  velocity,  potential  temperatures  and 
water  vapor  content  in  the  PBL,  n  is  a  unit  vector  normal  to  t,  h  is 
the  height  of  the  PBL,  Fu,  Fy,  FQ,  and  F^  are  functions  to  be  determined. 

Assuming  a  matching  layer  where  (1)  and  (2)  are  simultaneously 
valid,  one  can  obtain  a  relation  between  the  surface  fluxes  and  the 
average  PBL  parameters: 
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where  um  and  vm  are  the  components  of  the  meah  velocity  in  the  PBL  in 
t  and  n  directions,  respectively.  Universal  functions  A,  B,  C,  and  D 
determined  empirically.  It  is  generally  assumed  that  DwC  (see  Arya, 
1977). 

Using  equation  set  (3),  and  the  Wagara  observational  data, 
Yamada  (1976)  obtained  functions  A,  B,  and  C  (eqs.  (13)  through  (18) 
in  his  paper)  which  are  adopted  here  because  of  the  simple  functional 
form  and  accuracy. 

It  should  be  noted  that  in  Yamada 's  analysis  is  the  ver¬ 

tically  averaged  geostrophic  wind  in  the  PBL.  Here  Vm  is  the  verti¬ 
cally  averaged  mean  wind.  Thus  the  baroclinic  effects  are  minimized 
(Arya,  1978).  In  addition,  when  actually  applying  Eq.  (3)  in  a 
dynamic  model,  the  sum  of  the  squares  and  the  quotients  of  Eqs.  (3a) 
and  (3b)  are  used  to  evaluate  the  magnitude  and  direction  of  the 
surface  wind  stress. 
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III.  GENERALIZED  SIMILARITY  THEORY  AND  MULTI-LAYER  APPROACHES 

IN  A  ONE  DIMENSIONAL  TROPICAL  CYCLONE  MODEL 

As  stated  in  the  introduction,  the  parameterization  of  the  PBL 
based  on  the  GST  for  a  tropical  cyclone  model  is  very  appealing  because 
it  requires  only  one  layer  to  represent  the  PBL.  Thus  it  is  important 
to  test  the  universal  functions  decribed  by  Eq.  (3)  with  the  observa¬ 
tions  of  tropical  cyclones.  Unfortunately,  detailed  observations  of 
tropical  cyclone  PBL's  are  not  yet  available,  so  we  resort  to  use  a 
one-dimensional,  multi-layer  PBL  model  for  comparison.  Table  1  lists 
various  models  used  in  this  study. 

The  BCA  1-D  PBL  model  can  reproduce  observations  accurately  and 
it  has  been  used  successfully  in  parameterizing  the  PBL  in  numerical 
models  (Anthes  and  Chang,  1978;  Anthes  et  al . ,  1978).  The  model  has 
a  surface  layer  in  which  the  logarithmic-linear  relationships  for  wind 
and  temperature  (Businger  et  al.,  1971)  are  adopted.  It  should  be 
noted  that  the  same  logarithmic-linear  relationship  is  used  by  Yamada 
(1976)  to  generate  the  universal  functions.  Above  the  surface  layer, 
a  time  dependent  mixing  length  is  used.  The  mixing  length  is  assumed 
to  continuously  approach  its  asymptotic  value,  which  depends  on  the 
depth  and  stability  of  the  PBL.  The  e-folding  time  of  the  mixing 
length  is  proportional  to  the  size  of  the  energy-containing  eddies  in 
the  PBL  and  inversely  proportional  to  the  turbulent  kinetic  energy. 

The  ID  version  of  the  BCA  model  used  in  this  study  has  20  layers 
using  a  stretched  grid  from  the  surface  to  3  km.  The  lowest  layer 
defined  as  the  surface  layer,  has  a  thickness  of  25  m.  To  properly 
simulate  tropical  cyclone  condition,  a  rather  strong  pressure  gradient 
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TABLE  1 

Model  Used  For  Comparison 


Model s _ Type _ Dimensions _ PBL  Parameterization 


GST 

PBL 

1 

Single-layer  generalized 
similarity  theory 

BCA 

PBL 

1 

Multi-layer,  Busch  et  al. 

(1976) 

S 

Tropical  Cyclone 

2 

Single-layer,  generalized 
similarity 

D1 

Tropical  Cyclone 

2 

Single-layer,  Cn  =  0.0015 

CE  =  CH  =  0. 003u 

D2 

Tropical  Cyclone 

2 

Single-layer,  Cn  =  0.003 

CE  =  CH  =  0. 003u 

M 

Tropical  Cyclone 

2 

Multi-layer,  Busch  et  al . 
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force  is  specified  so  that  the  geostrophic  wind  speed  is  strongest  at  the 
surface  and  decreases  with  height  to  90%  of  its  surface  strength  at 
3  km.  The  initial  thermal  state  which  has  been  taken  from  the  model 
hurricane  of  Anthes  and  Chang  (1978)  shows  a  well  mixed  PBL  below  a 
capping  inversion  at  600  m  with  a  sharp  increase  of  potential  temperature 
and  decrease  of  specific  humidity  above  the  inversion.  The  initial  wind 
is  equal  to  the  geostrophic  wind  prescribed  for  each  run. 

The  BCA  model  was  integrated  12  hours  for  24  cases  with  the 
values  of  the  geostrophic  wind  set  at  15,  25  and  35  m  s”1  and  the  sea- 
surface  temperature  (SST)  varying  from  24.4  to  29.4°C.  These  variations 
produce  different  stability  conditions  and  large  rates  of  change  in  the 
initial  hours.  Diagnoses  using  Eq.  (3)  (or  GSA  model)  were  given  at 
model  times  of  1,  2,  3,  4,  6,  8,  10,  and  12  h.  The  vertically  inte¬ 
grated  u,  v,  0,  and  q  in  the  BCA  model  are  the  mean  PBL  properties  needed 
for  the  diagnoses.  Figure  1  shows  the  comparison  of  the  surface  wind 
stresses  between  the  diagnoses  from  GST  model  and  the  computation  of 
BCA  model.  The  GST  model  produced  nearly  the  same  results  (points  on 
the  diagonal  line  in  Figure  1)  as  the  BCA  model  for  small  u*  value  and 
for  the  unstable  conditions.  The  stress  diagnosed  by  the  GST  are 
higher  than  the  BCA  model  under  stable  conditions.  The  overestimate 
by  the  GST  model  lessens  as  instability  increases. 

The  larger  estimation  of  the  surface  stresses  for  stable  con¬ 
ditions  of  the  GST  model  is  due  to  the  greater  wind  shear  in  the  lower 
PBL  which  is  resolvable  in  the  BCA  model.  This  strong  wind  shear  is 
not  properly  represented  by  the  universal  functions  in  the  single  layer 
approach.  When  such  a  strongly  sheared  layer,  which  usually  accrmpanies 


a  low-level  inversion,  is  present,  the  PBL  mean  wind  may  not  be  a  proper 
velocity  scale  for  the  PBL.  Fortunately,  the  formation  of  such  a  low- 
level  inversion  is  not  likely  in  tropical  cyclones,  the  GST  model  there¬ 
fore  seems  to  estimate  the  momentum  fluxes  in  typical  tropical  cyclones 
quite  well. 

We  useda  constant  aerodynamic  drag  coefficient  of  1.5  x  10-3  to 
estimate  the  stresses  based  on  the  surface  winds  from  the  BCA  model. 
Because  the  drag  coefficient  is  independent  of  stability,  the  value  of 
the  stresses  depends  strongly  on  the  wind  speed.  A  comparison  with 
BCA  model  shows  that  most  of  the  values  cluster  along  the  three  hori¬ 
zontal  lines  and  scatter  within  the  vertical  bars  on  two  ends  of  each 
line  in  Figure  1.  These  scattering  and  clustering  points  are  also 
typical  for  the  drag  coefficient  model  in  heat  flux  and  moisture  flux. 

It  is  apparant  that  the  stresses  of  the  GSA  model  have  more  correlation 
and  less  scattering  to  those  of  the  more  sophisticated  BCA  model  as 
compared  to  the  aerodynamic  drag  coefficient  model. 

Figure  2  and  Figure  3  show  the  comparisons  of  the  upward  and 
downward  sensible  heat  fluxes,  respectively.  For  unstable  conditions. 
The  GST  model  underestimates  the  upward  heat  flux  because  of  the  for¬ 
mation  of  a  strong  superadiabatic  lapse  rate  region  near  the  surface 
in  the  BCA  model.  For  stable  conditions  the  GST  model  overestimates 
the  downard  heat  flux.  This  is  related  to  the  formation  of  low-level 
inversion  in  the  BCA  model  which  effectively  cuts  off  the  downward 
transport  of  heat.  These  characteristics  of  the  GST  model  are  found 
in  almost  all  single-layer  parameterizations  because  the  strong  gradient 
in  lower  PBL  cannot  be  resolved.  This  difficulty  cannot  be  easily 
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overcome  when  a  single- layer  PBL  has  to  be  employed  in  a  numerical  model 
for  economic  reasons. 

Figure  4  shows  the  comparison  of  the  moisture  flux  w 1  q '  •  The 
GST  model  gives  a  higher  w'q*  for  the  unstable  condition,  and  a  very 
similar  w'q'  for  near  neutral  conditions,  and  scattering  in  Fig.  4  is 
similar  to  that  in  Fig.  1,  the  more  symmetric  distribution  about  the 
diagonal  in  Fig.  4  is  presumably  related  to  the  larger  value  of 
in  the  lower  part  of  the  BCA  model. 

These  deviations  of  the  GST  model  from  the  BCA  model  are  much 
less  than  the  model  using  aerodynamic  drag  coefficients.  Using  a 
constant  drag  coefficient,  various  fluxes  depend  almost  solely  on  the 
wind  speed.  The  deviations  of  the  GST  model  from  the  BCA  model  is 
partially  inherited  from  the  scattering  of  the  original  observation 
data  from  which  the  universal  functions  are  deduced  (Yamada,  1976). 

The  comparison,  however,  does  indicate  a  distinct  behavior  of  the  GST 
model  for  different  stability.  It  should  be  pointed  out  that  the  GST 
has  been  applied  to  highly  non-stationary  conditions  in  this  analysis. 

We  note  that  as  the  PBL  approaches  a  quasi -steady  state,  the  GST  model 
diagnoses  becomes  comparable  to  the  BCA  model  results. 

The  experiments  with  the  one-dimensional  model  indicates  that 
the  GST  model  agrees  very  well  with  the  Busch  model  in  momentum  and 
moisture  flux  predictions.  The  GST  model  does  not  predict  accurate 
fluxes  when  there  are  strong  low-level  gradients  of  potential  temperature. 
The  surface  heat  flux,  however,  is  fortuitously  not  critically  important 
in  the  development  and  maintenance  of  tropical  cyclones  as  pointed  out 
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by  Rosenthal  (1971)  and  Anthes  and  Chang  (1978).  The  single  layer 
approach  with  GST  model  in  parameterizing  the  PBL  in  a  tropical  cyclone 
is  thus  adequate. 


IV.  INCORPORATION  WITH  AN  AXISYMMETRIC  TROPICAL  CYCLONE  MODEL 
A.  Model  Review 

The  ax i symmetric  tropical  cyclone  model  used  for  our  investi¬ 
gation  is  similar  to  the  one  described  in  Anthes  and  Chang  (1978) 
except  that  (1)  the  finest  resolution  at  small  radii  is  30  km,  and 
(2)  there  are  six  vertical  layers.  The  vertical  o-layers  are  bounded 
at  levels  o  =  0.0,  0.2,  0.3,  0.6,  0.8,  0.93,  and  1. 

Included  in  one  earlier  version  of  the  model  was  a  time  dependent 
equation  for  the  PBL  height  similar  to  that  in  Deardorff  (1972). 
Supplementary  experiments  showed  that  the  equation  predicts  infinite 
growth  of  PBL  near  the  center  of  the  cyclone  unless  large  horizontal 
diffusion  is  applied.  We,  therefore,  hold  the  height  of  PBL  constant 
between  a  =  0.93  and  1.  The  PBL  has  a  depth  of  650  =  700  m,  a  typical 
observed  PBL  depth  near  topical  disturbances  (Moss  and  Merceret,  1976; 
Nitta,  1974).  Charnock's  equation  (Delsol  et  al . ,  1970)  is  used  to 
compute  the  roughness  length.  The  initial  conditions  consist  of  a  weak 
vortex,  which  is  in  gradient  balance,  with  maximum  tangential  velocity  of 
17  m  s'1,  embedded  in  a  tropical  atmosphere  (Jordon,  1958).  The  lateral 
boundary  conditions  are  Dirichlet  for  the  thermodynamic  variables  and 
zero- vorti city  for  the  momentum. 

For  comparison  purposes  models  were  integrated  using  constant 
drag  coefficients  with  Cp  =  0.0015,  =  0.003  (model  Dl)  and  with 

Cp  =  0.003,  CH  *  CE  0.003  (model  D2).  A  nine-layer  tropical  cyclone 
model  (model  M),  identical  to  that  of  Anthes  and  Chang  (1978)  except 
for  the  30  km  horizontal  resolution,  is  also  integrated. 


B.  Development  into  Tropical  Cyclones 


The  development  from  the  initial  vortex  into  a  tropical  cyclone 
is  depicted  by  the  maximum  surface  wind  (Fig.  5)  and  minimum  pressure 
(Fig.  6).  After  the  initial  dissipation  due  to  friction.  Model  S 
rapidly  develops  into  an  intense  cyclone  and  reaches  a  quasi-steady 
state  at  about  40  hr.  The  maximum  surface  wind  is  about  57  m  s”1  and 
the  minimum  surface  pressure  is  approximately  958  mb. 

Models  D1  and  D2  also  intensify  after  the  initial  dissipation 
stage.  These  two  models  deviate  considerably  from  each  other  and  from 
Model  S  after  20  h.  Model  D1  reaches  a  final  intensity  of  35  m  s"1  and 
986  mb,  whereas  Model  D2  reaches  a  final  intensity  of  43  m  s-1  and  980 
mb.  Model  M  behaves  in  a  very  similar  fashion  as  Model  S, examplified 
by  the  continuing  intensification  between  20  to  40  h  and  the  final 
intensity  of  52  m  s"1  wind  and  956  mb  pressure. 

C.  Cyclone  Structure  at  Quasi -Steady  State 

The  tangential  velocity  field  at  48  h  for  Model  S  shows  a  strong 
and  concentrated  cyclonic  circulation  near  the  center  (Fig.  7).  The 
maximum  velocity  is  at  r  =  30  km  in  the  lowest  level  and  the  40  m  s”1 
contour  extends  upward  to  o  =  0.6.  The  tangential  circulation  diminishes 
upward  due  to  the  baroclinic  effect  of  the  warm  core  and  anti cyclonic 
circulation  occurs  at  the  o  =  0.2  level  outside  of  300  km.  The  inflow 
is  confined  to  the  PBL  with  maximum  radial  velocity  of  35  ms"1  just 
outside  the  region  of  the  maximum  tangetial  velocity  (Fig.  8).  Strong 
outflow  with  velocity  more  than  20  m  s"1  occurs  at  the  o  =  0.2  level. 
Except  for  the  shallow  inflow  and  outflow  layers,  the  radial  velocity 
is  small  in  mid-troposhere. 
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•(  Figure  9  shows  the  relative  humidity  (RH)  at  43  h  of  Model  S. 

The  RH  field  features  (1)  very  moist  PBL  and  outflow  canopy,  (2) 
relative  moist  convective  eyewall  region,  (3)  a  dry  eye  region  due  to 
j  the  strong  sinking  motion,  and  (4)  a  dry  mid-troposphere  outside  the 

eyewall  due  to  the  general  subsidence. 

The  temperature  field  of  the  quasi-steady  tropical  cyclone 
f  shows  a  warm  core  with  an  anomaly  of  more  than  +12°C.  The  +8°C  anomaly 

extends  to  r  =  300  km  in  the  outflow  layer.  The  +2°C  anomaly  contour 
between  a  =  0.06  and  a  =  0.9  outlines  the  subsidence.  The  cooler  PBL 

*  near  the  center  is  due  to  the  adiabatic  expansion  (and  therefore  cooling) 
in  the  inflow. 

0.  Air-Sea  Energy  Transfer 

*  We  will  use  budget  calculations  of  kinetic  energy,  moisture, 
and  heat  to  compare  model  results  and  reveal  the  importance  of  various 
physical  processes.  As  shown  in  Anthes  and  Chang  (1978),  the  budgets 
of  kinetic  energy  (K),  moisture  (Q)  and  enthalpy  (H)  for  a  tropical 
cyclone  model  can  be  expressed  schematically  by: 

*  f-"b  +  C  +  D¥  +  DL  +  DH 

H»Qb  +  EtP  +  QH 

) 

|f«Kb  +  A*Qp*Hs  *Hh 

where  the  terms  with  subscript  b's  denote  the  lateral  boundary  fluxes, 
and  terms  with  subscript  H's  denote  the  effects  due  to  horizontal 
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diffusion.  In  addition,  Dy,  and  D ^  represent  the  rates  of  kinetic 
energy  dissipation  due  to  surface  friction  and  cumulus  friction  re¬ 
spectively.  E  and  P  denote  evaporation  and  precipitation  rates 
respectively.  Qp  represents  latent  heat  release,  and  Hg  the  surface 
sensible  heat  transfer.  A  denotes  the  adiabatic  heating  rate  due  to 
vertical  motion;  C  denotes  the  kinetic  energy  conversion  rate  due  to 
cross-isobaric  flow. 

Within  the  model  domain  of  o  =  0  to  1  and  r  =  0  to  300  km, 

the  most  direct  contribution  of  the  PBL  processes  are  the  air-sea 

sensible  and  latent  heat  exchanges.  Figure  11  and  12  show  the  time 

series  of  Hs  and  E.  During  the  early  stages,  the  sensible  and  latent 

heat  transfers  of  D1  are  equal  or  greater  than  those  of  S.  The 

sensible  and  latent  heat  exchanges  in  S,  however,  grow  steadily  after 

12  -1 

10  hr  and  finally  reach  15  x  10  W  and  3  cm  d  respectively,  whereas 

in  Dl,  the  sensible  and  latent  heat  exchanges  only  grow  to  maxima  of 
12  -1 

7.5  x  10  W  and  1.8  cm  d  after  30  h.  After  the  maxima  are  reached, 
they  decrease  with  time  because  of  stabilization  and  saturation  of  the 
PBL. 

Apparently,  the  sensible  and  latent  heat  exchanges  in  S  are 
maintained  at  high  levels  primarily  due  to  their  nonlinear  dependence 
on  wind  speed  and  stability.  As  wind  speed  in  S  approaches  its  maximum 
after  30  h,  both  H$  and  E  are  enhanced  to  a  greater  extent  than  in  Dl. 
In  addition,  the  radial  circulation  in  S  is  also  increased  so  that  a 
slightly  unstable  and  dry  PBL  is  maintained  by  the  strong  inflow. 
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The  time  variation  of  E  for  model  Mis  very  close  to  model  S. 

The  sensible  heating  rate  H$,  however,  are  very  different.  This  is 
because  in  model  M  a  downward  heat  flux  across  the  PBL  top  is  allowed. 
Anthes  and  Chang  (1978)  argued  this  as  an  alternative  heat  source  for 
maintaining  the  quasi-isothermal  hurricane  PBL.  The  GSA  does  not  consider 
the  across  inversion  mising,  therefore  Model  S  depends  on  extracting  sensi¬ 
ble  heat  from  the  ocean  to  maintain  its  quasi-isothermal  inflow  layer. 

The  increased  radial  circulation  is  illustrated  by  the  con¬ 
version  rate  C  from  the  available  potential  energy  to  kinetic  energy 
shown  in  Fig.  13.  The  conversion  rate  is  approximately  proportional 
to  the  combined  strength  of  the  radial  circulation  and  the  warm  core. 

It  is  apparent  that  the  conversion  rates  in  Model  S  and  M  grow  rapidly 
between  10  to  20  h,  which  coincides  with  the  increased  evaporation. 

The  proportionally  faster  increase  of  C  than  H$  and  E  in  Model  S 
shows  that  the  increased  heat  transfer  creates  a  stronger  inflow, 
which,  in  turn,  maintains  a  stronger  heat  transfer  and  a  stronger 
inward  transport  of  the  dry,  cooler  ambient  air. 

The  two  different  PBL  parameterizations  also  gives  different 
dissipation  rates  (Fig.  14).  Note  that  Dy  in  S  is  smaller  before  15  h 
than  in  Dl,  and  it  eventually  asymptotes  to  three  times  stronger. 

Again,  this  is  due  to  the  nonlinear  dependence  of  the  surface  stress 
on  wind  speed. 

One  can  argue  that  the  nonlinear  dependence  of  the  air-sea 
energy  exchange  could  be  achieved  by  employing  wind-speed-dependent 
drag  coefficients.  However,  the  parameterization  based  on  similarity 
theory  is  superior.  The  heat  and  momentum  transfers  depend  not  only 


on  wind  speed  but  also  on  PBL  stability,  and  it  is  through  PBL  stability 
that  various  fluxes  are  correlated.  The  surface  roughness  appears  in 
the  formulation  as  a  length  scale  and  thus  the  parameterization  captures 
an  important  feedback  from  the  sea  state.  The  turning  of  wind  in  low 
levels,  neglected  in  the  parameterization  with  drag  coefficient,  it  also 
included  implicitly.  Therefore,  the  surface  stress  backs  from  the 
direction  of  the  PBL  mean  wind,  as  expected. 

Overall,  the  generalized  similarity  theory  is  quite  similar  in 
estimating  surface  fluxes  to  the  more  sophisticated  Busch  et  al  model  and 
appears  to  be  applicable  in  parameterization  of  the  PBL  in  the  tropical 
cyclone  models.  Such  parameterization  contains  several  realistic  PBL 
characteristics  such  as  the  mutual  dependence  of  various  fluxes  and  the 
backing  of  the  low-level  wind.  Although  it  cannot  resolve  the  strong 
gradient  in  a  stable  PBL  nor  the  across-inversion  mixing  in  some  single¬ 
layer  parameterizations,  it  is  still  superior  than  the  paremeterization 
using  stability-independent  drag  coefficient.  In  cases  that  a  multi-layer 
parameterization  is  economically  infeasible,  the  parameterization  outlined 
by  (3)  is  a  good  alternative. 
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V. 


SUMMARY 


A  PBL  parameterization  based  on  the  generalized  similarity 
theory  has  been  tested  for  tropical  cyclone  models.  We  conclude  that 
a  single  layer  approach  in  parameterizing  the  PBL  is  viable  although 
it  cannot  resolve  the  detailed  structure  of  the  PBL. 

This  parameterization  was  compared  with  the  ID  multi-level 
PBL  model  of  Busch  et  al.  (1976)  under  various  stabilities  and  highly 
nonstationary  conditions.  For  stationary  conditions,  the  computed  water 
vapor  flux,  agrees  well  with  the  Busch  model.  The  momentum  flux  is 
overestimated  and  the  sensible  heat  flux  is  underestimated  especially 
when  there  is  a  low-level  inversion. 

An  axi symmetric  tropical  cyclone  model  which  incorporated  this 
single-layer  parameterization  produced  very  realistic  dynamic  and 
thermodynamic  structures.  In  comparison  with  otherwise  identical 
tropical  cyclone  models  that  employ  constant  drag  coefficients,  the 
superiority  of  stress-dependent  drag  coefficients  in  the  development 
of  tropical  cyclones  was  demonstrated. 


* 
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FIGURES 


2  1 

1.  Surface  stresses  (m  s  )  computed  by  generalized  similarity  theory 
versus  those  computed  by  the  ID  PBL  model  for  stable  (•),  unstable  (o), 
and  neutral  (m)  conditions.  Points  should  fall  on  the  diagonal  line 

if  two  models  agree. 

2.  Same  as  Figure  1  except  for  upward  surface  heat  flux  in  unstable 
condition  (C  m  s-*). 

3.  Same  as  Figure  1  except  for  downward  surface  heat  flux  in  stable 
conditions  (C  m  s-i). 

4.  Same  as  Figure  1  except  for  surface  vapor  flux  (g  g“*  m  s-1). 

5.  Time  series  of  the  maximum  winds  for  models  S,  Dl,  D2,  and  M. 

6.  Time  series  of  the  minimum  surface  pressures  for  models  S,  Dl, 

D2,  and  M. 

7.  The  quasi-steady  tangential  velocity  (m  s"1)  of  model  S  at  48  h. 

Positive  values  denote  cyclonic  circulation,  negative  values  denote 
anti  cyclonic  circulation. 

8.  The  quasi-steady  radial  velocity  (m  s"*)  of  model  S  at  48  h. 

Positive  values  denote  outflow,  negative  values  denote  inflow. 

9.  The  quasi-steady  relative  humidity  [%)  of  Model  S  at  48  h. 

10.  The  quasi -steady  temperature  anomalies  (C)  of  Model  S  at  48  h. 

The  anomalies  are  defined  as  temperature  departures  from  the  initial 
state. 

11.  Time  series  of  surface  heat  flux  from  r  =  0  to  300  km  of  models  S, 

Dl,  and  M. 

12.  Same  as  Figure  11  except  for  averaged  evaporation  rate. 

13.  Same  as  Figure  11  except  for  KE  conversion  rate. 

14.  Same  as  Figure  11  except  for  KE  dissipation  rate  due  to  surface  friction. 


22 


(Busch  et  al)  X  10' 


Figure  1 


(Yamada)  X  10 


Central  Pressure  (mb) 


Figure 


Fipure 


Figure 


tw„uu*n  '•pu“'3  (w„on*a  ,w*‘01l° 


